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Abstract 

We have measured the differential production cross sections as a function of scaled mo- 
mentum x p = 2p/E cm of the identified hadron species 7r + , K + , K°, K*°, (ft, p, A , and of 
the corresponding antihadron species in inclusive hadronic Z° decays, as well as separately 
for Z° decays into light (u, d, s), c and b flavors. Clear flavor dependences are observed, 
consistent with expectations based upon previously measured production and decay proper- 
ties of heavy hadrons. These results were used to test the QCD predictions of Gribov and 
Lipatov, the predictions of QCD in the Modified Leading Logarithm Approximation with the 
ansatz of Local Parton-Hadron Duality, and the predictions of three fragmentation models. 
Ratios of production of different hadron species were also measured as a function of x p and 
were used to study the suppression of strange meson, strange and non-strange baryon, and 
vector meson production in the jet fragmentation process. The light-flavor results provide 
improved tests of the above predictions, as they remove the contribution of heavy hadron 
production and decay from that of the rest of the fragmentation process. In addition we 
have compared hadron and antihadron production as a function of x p in light quark (as 
opposed to antiquark) jets. Differences are observed at high x p , providing direct evidence 
that higher-momentum hadrons are more likely to contain a primary quark or antiquark. 
The differences for pseudoscalar and vector kaons provide new measurements of strangeness 
suppression for high-x p fragmentation products. 
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1 Introduction 



The production of jets of hadrons from hard partons produced in high energy collisions is 
believed to proceed in three stages. Considering the process e + e~ — > qq, the first stage 
involves the radiation of gluons from the primary quark and antiquark, which in turn may 
radiate gluons or split into qq pairs until their virtuality approaches the hadron mass scale. 
This process is in principle calculable in perturbative QCD, and three approaches have been 
taken so far: i) differential cross sections have been calculated [Q] for the production of up to 
4 partons to second order in the strong coupling a s , and leading order calculations have been 
performed recently for as many as 6 partons (see e.g. 0); ii) certain parton distributions 
have been calculated to all orders in a s in the Modified Leading Logarithm Approximation 
(MLLA) H; iii) "parton shower" calculations have been implemented numerically; these 
consist of an arbitrary number of q^-qg, g^gg and g—>qq branchings, with each branching 
probability determined from QCD in the Leading Logarithm Approximation. 

In the second stage these partons transform into "primary" hadrons. This "fragmenta- 
tion" process is not understood quantitatively and there are few theoretical predictions that 
do not explicitly involve heavy (c or b) quarks. Using perturbative QCD, Gribov and Lipatov 
have studied |5| the fragmentation of quarks produced in e + e~ collisions in the limit of high 
hadron momentum fraction x p = Phadron/Eb eam , and have related it to the proton structure 
function at high x = E quark / E proton . They predict that as x p — > 1 the distribution of x p 
for baryons is proportional to (1 — x p ) 3 , and that for mesons is proportional to (1 — x p ) 2 . 
Another approach is to make the ansatz of local parton-hadron duality (LPHD) [|J, that 
inclusive distributions of primary hadrons are the same, up to a normalization factor, as 
those for partons. Calculations using MLLA QCD, cut off at a virtual parton mass compa- 
rable with the mass of the hadron in question, have been used in combination with LPHD 
to predict that the shape of the distribution of £ = ln(l/x p ) for a given primary hadron 
species is approximately Gaussian within about one unit of the peak, that the shape can be 
approximated over a wider £ range by a Gaussian with the addition of small distortion terms, 
and that the peak position depends inversely on the hadron mass and logarithmically on the 
center-of-mass (cm.) energy. It is desirable to test the existing calculations experimentally 
and to encourage deeper theoretical understanding of the fragmentation process. 

In the third stage unstable primary hadrons decay into the stable particles that traverse 
particle detectors. This stage is understood inasmuch as proper lifetimes and decay branching 
ratios have been measured for many hadron species. However, these decays complicate 
fundamental fragmentation measurements because a sizable fraction of the stable particles 
are decay products rather than primary hadrons, and it is typically not possible to determine 
the origin of each detected hadron. Previous measurements at e + e~ colliders (see e.g. || |7j) 
indicate that decays of vector mesons, strange baryons and decuplet baryons produce roughly 
two-thirds of the stable particles; scalar mesons, tensor mesons and radially excited baryons 
have also been observed [0] , and there are large uncertainties on their contributions. Ideally 
one would measure every possible hadron species and distinguish primary hadrons from decay 
products on a statistical basis. A body of knowledge could be assembled by reconstructing 
heavier and heavier states, and subtracting their known decay products from the measured 
differential cross sections of lighter hadrons. 

Additional complications arise in jets initiated by heavy quarks, since the leading heavy 
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hadrons carry a large fraction of the beam energy, restricting that available to other pri- 
mary hadrons, and their decays produce a sizable fraction of the stable particles in the jet. 
Although decays of some B and D hadrons have been studied inclusively, there are large 
uncertainties in heavy hadron production, B° s and heavy baryon decay, and the suppression 
of gluon radiation from heavy quarks. The removal of heavy flavor events will therefore 
simplify the study of the fragmentation of light quarks into hadrons. 

A particularly interesting aspect of fragmentation is the question of what happens to the 
quark or antiquark that initiated the jet. A common prejudice is that the initial quark is 
"contained" as a valence constituent of a particular hadron, and that this "leading" hadron 
has on average a higher momentum than the other hadrons in the jet. The highly polarized 
electron beam delivered by the SLAC Linear Collider (SLC) gives a unique, high purity, 
unbiased tag of quark vs. antiquark jets, via the large electroweak forward-backward quark 
production asymmetry at the Z° resonance. We have previously observed || evidence for 
the production of leading baryons, K and K*°/K*° in light-flavor jets. The quantification 
of leading particle effects could lead to methods for identifying jets of specific light flavors, 
which could have a number of applications in ep and hadron-hadron collisions as well as in 
e + e~ annihilations. 

There are several phenomenological models of jet fragmentation, which combine mod- 
elling of all three stages of particle production; it is important to test their predictions. To 
simulate the parton production stage, the HERWIG @, JETSET and UCLA |I| event 
generators use a combination of first order matrix elements and a parton shower. To simulate 
the fragmentation stage, the HERWIG model splits the gluons produced in the first stage 
into qq pairs, and these quarks and antiquarks are paired up locally to form colorless clus- 
ters that decay into the primary hadrons. The JETSET model takes a different approach, 
representing the color field between the partons by a semi-classical string, which is broken, 
according to an iterative algorithm, into several pieces that correspond to primary hadrons. 
In the UCLA model, whole events are generated according to weights derived from the phase 
space available to their final states and the relevant Clebsch-Gordan coefficients. Each of 
these models contains arbitrary parameters that control various aspects of fragmentation and 
have been tuned to reproduce data from e + e~ annihilations. The JETSET model includes a 
large number of parameters that control, on average, the species of primary hadron produced 
at each string break, giving it the potential to model the observed properties of identified 
hadron species in great detail. In the HERWIG model, clusters are decayed into pairs of 
primary hadrons according to phase space, and the relative production of different hadrons 
is effectively governed by two parameters controlling the distribution of cluster masses. In 
the UCLA model, there is only one such free parameter, which controls the degree of locality 
of baryon-antibaryon pair formation. 

In this paper we present an analysis of tt^, K^, K°/K°, K*°/K*°, <f), p/p, and A°/A° 
production in hadronic Z° decays collected by the SLC Large Detector (SLD). The analy- 
sis is based upon the approximately 150,000 hadronic events obtained in runs of the SLC 
between 1993 and 1995. We measure differential production cross sections for these seven 
hadron species in an inclusive sample of hadronic Z° decays and use the results to test the 
QCD predictions of Gribov and Lipatov, the predictions of MLLA QCD+LPHD, and the 
predictions of the three fragmentation models just described, as well as to study the sup- 
pression of strange hadrons, baryons, and vector mesons in the fragmentation process. We 
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also measure these differential cross sections separately in Z° decays into light flavors (uu, 
dd and ss), cc and bb, which provide improved tests of the QCD predictions, new tests of 
the fragmentation models that separate the heavy hadron production and decay modelling 
from that of the rest of the fragmentation process, and cleaner measurements of strangeness, 
baryon and vector-meson suppression. In addition we update our measurements of hadron 
and antihadron differential cross sections in light quark jets, and use the results to make 
additional new tests of the fragmentation models and to make two new measurements of 
strangeness suppression at high x p . 

In section 2 we describe the SLD, including a detailed description of the Cherenkov Ring 
Imaging Detector, which is used to identify charged hadrons. In section 3 we describe the 
selection of hadronic events of different primary flavor, using impact parameters of charged 
tracks measured in the Vertex Detector, and the selection of light quark and antiquark hemi- 
spheres, using the large production asymmetry in polar angle induced by the polarization 
of the SLC electron beam. In section 4 we describe the hadron identification analyses and 
present results for flavor-inclusive events. In section 5 we present results separately for light- 
(Z° — ► uu,dd, ss), c- (Z° — > cc) and 6-flavor (Z° — > bb) events. In section 6 we use the 
flavor-inclusive and light-flavor results to test the QCD predictions of Gribov and Lipatov, 
and of MLLA QCD+LPHD. In section 7 we extract total production cross sections of each 
hadron species per hadronic event. In section 8 we update our measurements of leading 
particle production in light-flavor jets. In section 9 we present ratios of production of pairs 
of hadrons, and discuss the suppression of strange hadrons, baryons, and vector mesons in 
the fragmentation process. 



2 The SLD 



This analysis of data from the SLD |L2] used charged tracks measured in the Central Drift 
Chamber (CDC) [JD3] and silicon Vertex Detector (VXD) [TJJ], and identified in the Cherenkov 
Ring Imaging Detector (CRID) |15| . The CDC consists of 80 layers of sense wires arranged in 
10 axial or stereo superlayers between 24 and 96 cm from the beam axis. The outermost layer 
covers the solid angle range | cos#| < 0.68. The average spatial resolution for hits attached to 
charged tracks is 92 fim. Momentum measurement is provided by a uniform axial magnetic 
field of 0.6 T. The VXD and CRID are described in the following subsections. 

Energy deposits reconstructed in the Liquid Argon Calorimeter (LAC) [|16) were used in 



the initial hadronic event selection and in the calculation of the event thrust []T7| axis. The 
LAC is a lead-liquid argon sampling calorimeter covering the solid angle range | cos#| < 0.98, 
which is segmented into 33x36 mrad projective towers, each comprising two electromagnetic 
sections and two hadronic sections, for a total thickness of 2.8 interaction lengths. The 
energy resolution is measured to be a = 15%VE for electromagnetic showers and 60%\/E 
for hadronic showers, where E is the energy in GeV. 

2.1 The SLD Vertex Detector 

Flavor tagging of events for this analysis was accomplished with the original SLD Vertex 



Detector [fn}| ; which was composed of 480 charge-coupled devices containing a total of 120 
million 22x22 /im 2 pixels, arranged in four concentric layers of radius between 2.9 and 4.2 
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Figure 1: Distribution of transverse impact parameters of tracks in e + e~ 
with respect to the primary interaction point measured in hadronic events. 
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cm. The outermost layer covered the solid angle range |cos#| < 0.75, and the azimuthal 
arrangement was such that a track would always encounter one of the two innermost layers 
and one of the two outermost layers; the average number of reconstructed hits per track was 
2.3. The 3-D spatial resolution for these hits was measured to be 5.5 ^m. 

Here we used only the information in the plane transverse to the beam axis. The im- 



pact parameter resolution in this plane was measured |L8| from the distribution of miss 
distances between the two tracks in Z° — > /U + /i~ events to be 11 /zm for 45.6 GeV/c muons 
reconstructed including at least one hit in the VXD. The transverse position of the pri- 
mary interaction point (IP) was measured using tracks in sets of ~30 sequential hadronic 
Z° decays, with a resolution measured from the distribution of impact parameters in the 
statistically independent /i-pair event sample (see fig. |l|) of 7 ± 2 /im. The impact parameter 
resolution for lower momentum tracks was determined using tracks in hadronic Z° decays, 
corrected for the contributions from decays of heavy hadrons. Including the uncertainty on 
the IP, a resolution of ll©70/(p_i_ sin 3 / 2 9) /jm was obtained, where p± is the track momen- 
tum transverse to the beam axis in GeV/c and 9 is the polar angle of the track with respect 
to the beam axis. 



2.2 The SLD Cherenkov Ring Imaging Detector 

Identification of charged tracks is accomplished with the barrel CRID ]n| , which covers the 
solid angle range | cos 6"| < 0.68. Through the combined use of liquid C 6 F 14 and gaseous 
C5F12+N2 radiators, the barrel CRID is designed to perform efficient separation of charged 
pions, kaons and protons over most of the momentum range in e + e~ annihilations at the 
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Z°, 0.3 < p < 46 GeV/c. A charged particle that passes through a radiator of refractive 
index n with velocity (3 above Cherenkov threshold, (3 > (3q = 1/n, emits photons at an 
angle 6 C = cos" 1 {1 / (3n) with respect to its flight direction. In the SLD, a charged particle 
exiting the CDC encounters a 1 cm thick liquid radiator, contained in one of 40 radiator 
trays. If the momentum of the particle is above its liquid Cherenkov threshold, UV photons 
are emitted in a cone about the particle flight direction. This 1-cm thick cone expands over 
a standoff distance of ~12 cm and each photon can enter one of 40 time projection chambers 
(TPCs) through an inner quartz window. 

The TPCs contain a photosensitive gas, ethane with ~0.1% TMAE [[0|. The resulting 
single photoelectrons drift along the beam direction to a wire chamber where the conversion 
point of each Cherenkov photon is measured in three dimensions using drift time, wire 
address and charge division. These positions are used to reconstruct a Cherenkov angle with 
respect to the extrapolated charged track. Liquid rings span 2-3 TPCs in azimuth and can 
be split between TPCs in the forward and backward hemispheres. 

The particle may then continue through a TPC, where it ionizes the drift gas, saturating 
the readout electronics, which were designed for single-electron detection, on 2-7 anode wires 
and effectively deadening ~5 cm 2 of detection area. Following the TPC, the particle passes 
through ~40 cm of the gas radiator volume. Radiated Cherenkov photons are focussed by 
one of 400 spherical mirrors onto the outer quartz window of a TPC. Gas rings are typically 
2.5 cm in radius at the TPC surface, and the mirrors are positioned such that no ring is 
focussed near an edge of a TPC or near the region saturated by its own track. The mirror 
arrangement and the large size of the liquid rings make the identification performance largely 
independent of the proximity of the track to any jet axis. 

The average liquid (gas) Cherenkov angle resolution was measured from the data to be 16 
(4.5) mrad, including the effects of residual misalignments of the TPCs, radiator trays and 
mirrors, and track extrapolation resolution. The local or intrinsic resolution was measured 
to be 13 (3.8) mrad, consistent with the design value. The average number of detected 
photons per full ring for tracks with (3 — 1 was measured in /i-pair events to be 16.1 (10.0). 
For hadronic events, a set of cuts was applied to reduce backgrounds from spurious hits 
and cross-talk from saturating hits, resulting in an average of 12.8 (9.2) accepted hits per 
ring. The average reconstructed Cherenkov angle for (3 = 1 tracks was 675 (58.6) mrad, 
corresponding to an index of refraction of 1.281 (1.00172), and Cherenkov thresholds of 0.17 
(2.4) GeV/c for charged pions, 0.62 (8.4) GeV/c for kaons and 1.17 (16.0) GeV/c for protons. 
This index was found to be independent of position within the CRID and the liquid index 
was found to be constant in time. Time variations in the gas index of up to ±0.00007 were 
tracked with an online monitor and verified in the data. 

Tracks were identified using a likelihood technique |19j . For each of the five stable charged 
particle hypotheses i = e, p, n, K, p, a likelihood Lj was calculated based upon the number 
of detected photoelectrons and their measured angles, the expected number of photons, the 
expected Cherenkov angle, and a background term. The background included the effects 
of overlapping Cherenkov radiation from other tracks in the event as well as a constant 
term normalized to the number of hits in the TPC in question that were not associated 
with any track. Particle separation was based upon differences between logarithms of these 
likelihoods, Li = In Li. 

The particle identification performance of the CRID depends on the track selection and 
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likelihood difference requirements for a given analysis. Here we discuss the example of the 
hadron fractions analysis described in section 4.1, where we consider only the three charged 
hadron hypotheses i = 7i,K,p. For tracks with p < 2.5 (p > 2.5) GeV/c, a particle was 
identified as species j if Cj exceeded both of the other log-likelihoods by at least 5 (3) units. 
We quantify the performance in terms of a momentum-dependent identification efficiency 
matrix E, each element of which represents the probability that a selected track from a 
true i-hadron is identified as a j-hadron, with i, j = ir,K,p. The elements of this matrix were 
determined where possible from the data p0 |. For example, tracks from selected K® and 



t decays were used as "pion" test samples, having estimated kaon plus proton contents of 
0.3% and 1.7% respectively. Figure |2| shows the probability for these tracks to be identified 
as pions, kaons and protons as a function of momentum. Also shown are results of the 
same analysis of corresponding samples from a detailed Monte Carlo (MC) simulation of the 
detector. The MC describes the momentum dependence well and reproduces the efficiencies 
to within ±0.03. Functional forms were fitted to the data, chosen to describe the momentum 
dependence of both data and simulated test samples, as well as that of simulated true pions 
in hadronic events. The simulation was used to correct the fitted parameters for non-pion 
content in the K® and r samples and differences in tracking performance between tracks 
in these samples and those from the IP in hadronic events. The resulting identification 
efficiency functions, E nn , E n K and E np , are shown in the leftmost column of fig. [3]. 

A similar procedure using only tt and p likelihoods was used to measure the 7r-p separation 
in the liquid (gas) system for p > 2 (17) GeV/c, and the simulation was used to convert that 
into E pp , shown in the bottom right of fig. [3[ E pp over the remaining momentum range, 
as well as the tt-K separation in the gas system below and near kaon threshold (p < 10 
GeV/c), was measured using protons from decays of tagged lambda hyperons [p0| . The 
remaining efficiencies in fig. § were derived from those measured, using the simulation. For 
example, Ekk is equal to E nn for momenta in the ranges 1.5 < p < 2.5 and 15 < p < 25 
GeV/c, since both species are well above the relevant Cherenkov threshold and their expected 
Cherenkov angles differ from that of the proton by an amount large compared with the 
angular resolution. Outside these ranges, E KK was related to E nn by a function derived 
from the simulation to account for the effects of the reduced photon yield near the kaon 
Cherenkov threshold and the fact that the expected kaon ring radius lies between those of 
the pion and proton. 

The bands in fig. ^| encompass the upper and lower systematic error bounds on the 
efficiencies. The discontinuities correspond to the 7r and K Cherenkov thresholds in the gas 
radiator. For the diagonal elements, the systematic errors correspond to errors on the fitted 
parameters and are strongly positively correlated across each of the three momentum regions. 
For the off-diagonal elements, representing misidentification rates, a more conservative 25% 
relative error was assigned at all points to account for the limited experimental constraints 
on the momentum dependence. These errors are also strongly positively correlated among 
momenta. The identification efficiencies in fig. |3] peak near or above 0.9 and the pion 
coverage is continuous from 0.3 GeV/c up to approximately 35 GeV/c. There is a gap in 
the kaon-proton separation between about 7 and 10 GeV/c due to the limited resolution of 
the liquid system and the fact that neither species is far above Cherenkov threshold in the 
gas system. The proton coverage extends to the beam momentum. Misidentification rates 
are typically less than 0.03, with peak values of up to 0.07. 
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Figure 2: Efficiencies for selected tracks from K® (squares) and r (circles) decays to be 
identified as each hadron species in the CRID. The solid symbols represent the data and the 
open symbols the simulation. 
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Figure 3: Calibrated identification efficiencies for tracks used in the charged hadron fractions 
analysis. The half-widths of the grey bands represent the systematic uncertainties, which are 
strongly correlated between momenta. Note the expanded vertical scale for the off-diagonal 
efficiencies. 
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Efficienf 

uu, dd, ss 


;y for Z 

CC 




bb 


Con 

uu, dd, ss 


nositioT 

CC 


bb 


light-tag 

c-tag 

6-tag 


0.845 
0.153 
0.002 


0.438 
0.478 
0.084 


0.075 
0.331 
0.594 


0.849 
0.378 
0.009 


0.124 
0.333 
0.100 


0.027 
0.290 
0.891 



Table 1: Efficiencies for simulated events in the three flavor categories to be tagged as 
light, c or b events. The three rightmost columns indicate the composition of each simulated 
tagged sample assuming the Standard Model relative flavor production. 



3 Event Selection 



The trigger and initial selection of hadronic events are described in | 2~T| . The analysis 
presented here is based on charged tracks measured in the CDC and VXD. A set of cuts 
was applied in order to select events well-contained within the detector acceptance. Tracks 
were required to have (i) a closest approach to the beam axis within 5 cm, and within 10 cm 
along the beam axis of the measured IP, (ii) a polar angle 9 with respect to the beam axis 
with | cos#| < 0.80, (iii) a momentum transverse to this axis p± > 150 MeV/c, and (iv) a 
momentum p < 50 GeV/c. Events were required: to contain a minimum of seven such tracks; 
to contain a minimum visible energy E vis > 18 GeV, calculated from the accepted tracks, 
assigned the charged pion mass; to have a thrust axis polar angle 9 t with respect to the 
beam axis, calculated from calorimeter clusters, with | cos9 t \ < 0.71; and to have good VXD 
data |T!| and a well- measured IP position. A sample of 90,213 events passed these cuts. For 



the analyses using the CRID, the additional requirements were made that the CRID high 
voltage was on and that there was a good drift velocity measurement, resulting in a sample 
of 79,711 events. The non-hadronic background was estimated to be 0.1%, dominated by 
Z° — > t + t~ events. 

Samples of events enriched in light and b primary flavors were selected based on signed 
impact parameters 5 of charged tracks with respect to the IP in the plane transverse to the 
beam. For each event we define n S i g to be the number of tracks passing a set of impact- 



parameter quality cuts [|1J| that have impact parameter greater than three times its estimated 
error, 5 > 3as- Events with n sig = were assigned to the light-tagged sample and those with 
n sig > 3 were assigned to the 6-tagged sample. The remaining events were classified as a 
c-tagged sample. The light-, c- and 6-tagged samples comprised 60.4%, 24.5% and 15.2% of 
the selected hadronic events, respectively. The tagging efficiencies and sample purities were 
estimated from our Monte Carlo simulation and are listed in table p]. 

Separate samples of hemispheres enriched in light-quark and light-antiquark jets were 
selected from the light-tagged event sample by exploiting the large electroweak forward- 
backward production asymmetry with respect to the beam direction. The event thrust axis 
was used to approximate the initial qq axis and was signed such that its ^-component was 
along the electron beam direction, t z > 0. Events in the central region of the detector, where 
the production asymmetry is small, were removed by the requirement \t z \ > 0.2, leaving 74% 
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of the light-tagged events. The quark-tagged hemisphere in events with left- (right-)handed 
electron beam polarization was defined to comprise the set of tracks with positive (negative) 
momentum projection along the signed thrust axis. The remaining tracks in each event 
were defined to be in the antiquark-tagged hemisphere. For the selected event sample, the 
average magnitude of the polarization was 0.73. Using this value and assuming Standard 
Model couplings, a tree- level calculation gives a quark (antiquark) purity of 0.73 in the 
quark- ( ant iquark- ) t agged sample . 



4 Hadron Identification Analysis 

In the following subsections we discuss details of the analysis for three categories of identified 
hadrons: charged tracks identified as n ± , K ± or p/p in the CRID; K® and A°/A° recon- 
structed in their charged decay modes and tagged by their long flight distance; and K*°/K*° 
and <p reconstructed in charged decay modes including one and two identified K^, respec- 
tively. The resulting differential cross sections for these seven hadron species in inclusive 
hadronic Z° decays are presented in the last subsection. 



4.1 Charged Hadron Fractions 

Reconstructed charged tracks were identified as charged pions, kaons or protons using infor- 
mation from only the CRID liquid (gas) radiator for tracks with p < 2.5 (p > 7.5) GeV/c; 
in the overlap region, 2.5 < p < 7.5 GeV/c, liquid and gas information was combined. Ad- 



ditional track selection cuts |20| were applied to remove tracks that interacted or scattered 
through large angles before exiting the CRID and to ensure that the CRID performance was 
well-modelled by the simulation. Tracks were required to have at least 40 CDC hits, at least 
one of which was at a radius of at least 92 cm, to extrapolate through an active region of the 
appropriate radiator(s), and to have at least 80 (100)% of their expected liquid (gas) ring 
contained within a sensitive region of the CRID TPCs. The latter requirement included re- 
jection of tracks with p > 2.5 GeV/c for which there was a saturated CRID hit within a 5 cm 
radius (twice the maximum ring radius) of the expected gas ring center. Tracks with p < 7.5 
GeV/c were required to have a saturated hit within 1 cm of the extrapolated track, and 
tracks with p > 2.5 GeV/c were required to have either such a saturated hit or the presence 
of at least four hits consistent with a liquid ring. These cuts accepted 47%, 28% and 43% of 
the tracks within the CRID acceptance in the momentum ranges p < 2.5, 2.5 < p < 7.5 and 
p > 7.5 GeV/c, respectively. For momenta below 2 GeV/c, only negatively charged tracks 
were used in order to reduce the background from protons produced in particle interactions 
with the detector material. 

In each momentum bin we measured the fractions of the selected tracks that were identi- 
fied as pions, kaons and protons. The observed fractions were related to the true production 
fractions by an efficiency matrix, composed of the values shown in fig. || This matrix was 
inverted and used to unfold our observed identified hadron fractions. This analysis proce- 
dure does not require that the sum of the charged hadron fractions be unity; instead the 
sum was used as a consistency check, which was found to be satisfied at all momenta (see 
fig. |j). In some momentum regions we cannot distinguish two of the three hadron species, 
so the procedure was reduced to a 2x2 matrix analysis and we present only the fraction of 
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the identified species, i.e. protons above 35 GeV/c and pions below 0.75 GeV/c and between 
7.5 and 9.5 GeV/c. 

Electrons and muons were not distinguished from pions; this background was estimated 
from the simulation to be about 5% of the tracks in the inclusive flavor sample, predominantly 
from c- and 5-flavor events. The fractions were corrected using the simulation for the lepton 
backgrounds, as well as for the effects of beam-related backgrounds, particles interacting in 
the detector material, and particles decaying outside the tracking volume. The conventional 
definition of a final-state charged hadron was used, namely a charged pion, kaon or proton 
that is either from the primary interaction or a direct decay product of a hadron that has 
proper lifetime less than 3xl0 _10 s and is itself a primary or a decay product of a primary 
hadron. 

The measured charged hadron fractions in inclusive hadronic Z° decays are shown in 
fig. |] and listed in tables Q-|4]. The systematic errors were determined by propagating the 
errors on the calibrated efficiency matrix (see sec. 2.2) and correspond to uncertainties in the 
average number of photons detected per track and the average resolution on the measured 
Cherenkov angles. They are therefore strongly positively correlated across each of the three 
momentum regions, p < 2.5, 2.5 < p < 7.5 and p > 7.5 GeV/c, and are indicated by the 
pairs of dashed lines in fig. |j. The errors on the points below ~6 GeV/c are dominated by 
the systematic uncertainties; for the points above ~15 GeV/c the errors have roughly equal 
statistical and systematic contributions. 

Pions are seen to dominate the charged hadron production at low momentum, and to 
decline steadily in fraction as momentum increases. The kaon fraction rises steadily to about 
one-third at high momentum. The proton fraction rises to a plateau value of about one-tenth 
at about 10 GeV/c. Where the momentum coverage overlaps, these measured fractions were 
found to be consistent with an average of previous measurements at the Z° [22, |23|, [24| . 
Measurements based on ring imaging and those based on ionization energy loss rates cover 
complementary momentum ranges and can be combined to provide continuous coverage over 
the range 0.22 < p < 45.6 GeV/c. 

Differential production cross sections were obtained by multiplying these fractions by our 
measured inclusive charged particle differential cross section, corrected, using our simulation, 
for the contribution from leptons. The integral of this cross section was constrained to be 
20.95 tracks per event, an average [53] of charged multiplicity measurements in Z° decays, 
and the momentum-dependence of our track reconstruction efficiency was checked by com- 
paring the momentum distributions of charged tracks in data and simulated r ± decays. We 
include a 1.7% error on the average multiplicity as a systematic normalization uncertainty, 
as well as a momentum-dependent uncertainty of O.llx \p — 3.8 GeV/c|%, derived from the 
study of r ± decays. The inclusive charged particle differential cross section is listed in table 
[5], and the resulting differential cross sections per hadronic event per unit x p for the identified 
hadrons are listed in tables 0-41[ The 1.7% normalization uncertainty is not included in the 
systematic error listed for any of the identified hadrons, nor is it included in the error bars 
in any of the figures. 
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0.4 1 4 10 40 

4-98 

8400A6 Momentum (GeV/c) 

Figure 4: Measured charged hadron production fractions in hadronic Z° decays. The circles 
represent the 71^ fraction, the squares the K ± fraction, the triangles the p/p fraction, and 
the open circles the sum of the three fractions. The error bars in the upper plot are statistical 
only; the dashed lines indicate the systematic errors, which are strongly correlated between 
momenta. The error bars on the sum are statistical and systematic added in quadrature. 
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Xp Range 


<x p > 


u 


1/N dn n /dx p 


o 


008- 


-0 


010 


o 


009 


o 


963±0 


004±0 


014 


482 3±2 


3 


±7 


2 


o 


010- 


-0 


012 


o 


Oil 


o 


924±0 


004±0 


006 


439 0±2 


3 


±3 


7 


o 


012- 


-0 


014 


o 


013 


o 


921±0 


003±0 


006 


400 5±2 

x V ' V ' > t_ / l ^ 


o 


±3 


3 


o 


014- 


-0 


016 


o 


015 


o 


906±0 


004±0 


006 


356. 1±1 


9 


±3 


o 


o 


016- 


-0 


022 


o 


019 


o 


886±0 


002±0 


006 


292 8±1 

— 1...' 1 L 


o 


±2 


4 


o 


022- 


-0 


027 


o 


025 


o 


872±0 


003±0 


006 


228 5±1 


o 


±1 


9 





027- 


-0 


033 


o 


030 


o 


831±0 


003±0 


006 


176. 6±0 


9 


±1 


4 





033- 


-0 


038 


o 


036 


o 


820±0 


004±0 


006 


144.4±0 


8 


±1 


2 


o 


038- 


-0 


044 


o 


041 


o 


823±0 


004±0 


010 


121. 7±0 


8 


±1 


6 


o 


044- 


-0 


049 


o 


047 


o 


806±0 


006±0 


015 


102 5±0 


9 


±1 


g 


o 


049- 


-0 


055 


o 


052 


o 


812±0 


008±0 


020 


89.2±0 


9 


±2 


2 


o 


055- 


-0 


060 


o 


058 


o 


788±0 


007±0 


029 


75.3±0 


8 


±2 


8 


o 


060- 


-0 


066 


o 


063 


o 


779±0 


007±0 


016 


66.0±0 


7 


±1 


4 


o 


066- 


-0 


071 


o 


069 


o 


763±0 


007±0 


010 


57.81±0 


60 


±0 


81 





071- 


-0 


077 


o 


074 


o 


767±0 


007±0 


009 


51 63±0 


56 


±0 


60 


o 


077- 


-0 


082 


o 


079 


o 


761±0 


007±0 


009 


45 95±0 


52 


±0 


54 


o 


082- 


-0 


088 


o 


085 


o 


750±0 


007±0 


008 


41.35±0 


49 


±0 


49 





088- 


-0 


099 


o 


093 


o 


743±0 


006±0 


008 


35 24±0 


32 


±0 


42 





099- 


-0 


110 


o 


104 


o 


714±0 


006±0 


008 


28 12±0 


29 


±0 


35 


o 


110- 


-0 


121 


o 


115 


o 


705±0 


007±0 


009 


23 57±0 


27 


±0 


30 


o 


121- 


-0 


143 


o 


131 


o 


695±0 


005±0 


009 


18 32±0 


17 


±0 


24 


o 


143 


-0 


164 


o 


153 


o 


670±0 


006±0 


009 


13 22±0 


14 


±0 


19 


o 


164- 


-0 


186 


o 


175 


o 


651±0 


006±0 


009 


9.84±0 


11 


±0 


15 


o 


186- 


-0 


208 


o 


197 


o 


644±0 


007±0 


008 


7.47±0 


09 


±0 


11 


o 


208- 


-0 


230 


o 


219 


o 


625±0 


008±0 


007 


5.711±0 


083±0 


080 


o 


230- 


-0 


252 


o 


241 


o 


611±0 


009±0 


006 


4.414±0 


074±0 


063 


o 


252- 


-0 


274 


o 


263 


o 


618±0 


010±0 


010 


3 612±0 


06^ 


3±0 


072 


o 


274- 


-0 


296 


o 


285 


o 


608±0 


011±0 


010 


2 886±0 


061±0 


060 


o 


296- 


-0 


318 


o 


307 


o 


583±0 


012±0 


Oil 


2 206±0 


054±0 


049 


o 


318- 


-0 


351 


o 


334 


o 


578±0 


012±0 


012 


1.739±0 


040±0 


044 


o 


351- 


-0 


384 


o 


366 


o 


603±0 


014±0 


015 


1 350±0 


036±0 


040 





384- 


-0 


417 





400 





523±0 


017±0 


016 


0.874±0 


031±0 


032 





417- 


-0 


450 





432 





520±0 


021±0 


020 


0.670±0 


029±0 


029 





450- 


-0 


482 





465 





534±0 


024±0 


024 


0.520±0 


026±0 


025 





482- 


-0 


526 





503 





508±0 


028±0 


027 


0.355±0 


021±0 


020 





526- 


-0 


570 





547 





514±0 


036±0 


031 


0.248±0 


018±0 


016 





570- 


-0 


658 





609 





501±0 


040±0 


038 


0.146±0 


012±0 


012 





658- 


-0 


768 





704 





580±0 


076±0 


053 


0.071±0 


009±0 


007 



Total Observed/Evt. 



14.52±0.02 ±0.27 



Table 2: Charged pion fraction f n and differential cross section (l/N)dn 7r /dx p per hadronic 
Z° decay. <x p > is the average x p - value of charged tracks in each bin. The last row gives 
the integral over the x p range of the measurement. The first error is statistical, the second 
systematic. A 1.7% normalization uncertainty is included in the systematic error on the 
integral, but not in those on the cross section. 
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Xp Range 


<x p > 


fK 


1/N dn K /dx 


V 





016- 


-0 


022 


n 




n i n 

019 


n 




n/ 1 t i n 

067±0 


nn 1 i n 

001±0 


nno 

002 


oo oo i n 

22.28±0 


A T 

47 


i n 

±0 


£0 

oo 


n 




022- 


-0 


027 


n 




nor 

025 


n 




no 1 i n 

081±0 


nno i n 

002±0 


nno 

002 


oi oo i n 

21.22±0 


A r 

45 


i n 

±0 


^o 
02 


n 




A O ^7 

027- 


-0 


noo 

033 


n 




non 

030 


n 




nnn i n 

090±0 


nno i n 

002±0 


nno 

003 


i n i n i n 

19.10±0 


A O 

43 


i n 

±0 


04 


n 




066- 


A 




noo 

038 


n 




no/ 1 

036 


n 




i no i n 

102±0 


nno i n 

002±0 


nn r 

005 


i o no i n 

18.02±0 


A O 

43 


i n 

±0 


on 
80 





non 

068- 


A 

-0 


n a a 

044 


n 




n a 1 

041 


A 




111 in 

111±0 


nno i n 

003±0 


nn/ 1 

006 


1 z 1 A r in 

16.45±0 


A r 

45 


i n 

±0 


94 


n 




C\ A A 

044- 


A 

-0 


r\ Af\ 

049 


n 




n a t 

047 


A 




i ot i n 

127±0 


nn /i i n 

004±0 


nno 

008 


i z 1 i o i n 

16.13±0 


A n 

49 


i 1 

±i 


no 
Oo 


n 


049- 


A 

-0 


nrr 

055 


n 




n r o 

052 


A 




i o t i n 

127±0 


nn r i n 

005±0 


n i n 

010 


i o no i n 

13.98±0 


r o 

53 


i i 

±i 


14 


n 




nrr 

055- 


A 

-0 


n/^n 

060 


n 




n r o 

058 


A 




i o r in 

125±0 


nn/ 1 i n 

006±0 


noo 

022 


1 1 n/ 1 i n 

11.96±0 


r ^ 

54 


1 o 

±2 


11 


n 
U 


UOU- 


n 
-u 


UOO 


n 
U 


Uoo 


n 
U 


1 on i n 
10U±U 


nn<3_i_n 
UU0±U 


Ulo 


1 i no i n 
ll.U0±U 


A O 


i i 
±1 




n 
U 


Uoo- 


n 

-0 


0/1 


n 



n^n 
U09 


n 



1 rn i n 
150±0 


nn^_i_n 
UUo±U 


012 


ll.o/±U 


4o 


i n 
±0 


O ( 





0/1- 


n 

-0 


M 


n 



0(4 


n 



1 on i n 
139±0 


00 ( ±0 


mo 
012 


n oo i n 

9.38±0 


44 


±0 


70 


n 



077- 


-0 


082 


n 




n^7n 

079 





1 c ^7 i n 

157±0 


nn^7 i n 

007±0 


mo 

013 


n c i i n 

9.51±0 


44 


±0 


/ 


a 
(J 


Uoz- 


(J 


Uoo 


A 
(J 


nor 


A 
(J 


i r ^7 i n 

157±U 


nnoj_n 
(JUczbU 


n i o 

Ulo 


o /^o in 
0.00±U 


44 


±0 


/ 2 


U 


Uoo- 


-U 


nnn 

uyy 


n 

u 


nno 

uyd 


A 
U 


1 /? o i n 

loo±U 


UUf±U 


n 1 /i 
U14 


1-7 n/? i n 

( .yo±u 


31 


±0 


OO 


n 
U 


nnn 

uyy- 


-U 


i in 
111) 


n 

u 


1U4 


A 

u 


1 o7_l n 
lo/±U 


uuy±i) 


ni £ 
Ulo 


T 0^7 i n 


34 


±0 


00 


U 


i in 

11U- 


-U 


1 Ol 

Izl 


n 
U 


llo 


A 
U 


ono i n 

zUz±U 


m 1 i n 
U11±U 


Ulo 


T/l i n 
0. /4±U 


37 


±0 


fin 
ou 


a 

o 


Izl- 


A 

-0 


1 A O 

143 


A 




1 O 1 

131 


A 




1 nn i n 

iyy±u 


mi i n 

011±0 


noo 

023 


r o a i n 

o.24±U 


29 


±0 


01 


a 

u 


143- 


A 

-0 


1 /? /I 

lo4 


A 




i r o 

153 


A 




OA ^7 i n 

207±0 


non i n 

020±0 


n /i 1 
U41 


/I no in 

4.08±0 


40 


±0 


qa 
oU 


A 


o n o 

208- 


A 

-0 


o on 

230 


A 




o i n 

219 


A 




or/ 1 i n 

256±0 


nnn i n 

009±0 


noo 

033 


o o /i in 

2.34±0 


08 


±0 


on 
3U 


A 




230- 


A 

-0 


252 


A 




O A 1 

241 


A 




opn i n 

269±0 


nnn i n 

009±0 


007 


1 C\ AT \ C\ 

1.947±0 


065±0 


Oo / 


A 


o r o 

2o2- 


A 

-0 


274 


A 




263 


A 




ot /i i n 

274±0 


nnn i n 

009±0 


007 


1 /?no i n 

1.603±0 


057±0 


(\A O 

04z 


A 


274- 


A 

-0 


on/"* 

296 


A 




o o r 

285 


A 




o i n 

270±0 


n 1 n i n 

010±0 


nn/ 1 

006 


i o o i in 

1.281±0 


050±0 


no /i 
0o4 





296- 


-0 


318 





307 





298±0 


011±0 


007 


1.127±0 


045±0 


non 
OoO 





318- 


-0 


351 





334 





310±0 


011±0 


008 


0.933±0 


034±0 


027 





351- 


-0 


384 





366 





299±0 


012±0 


009 


0.669±0 


029±0 


023 





384- 


-0 


417 





400 





324±0 


015±0 


012 


0.541±0 


026±0 


023 





417- 


-0 


450 





432 





383±0 


019±0 


016 


0.493±0 


026±0 


023 





450- 


-0 


482 





465 





366±0 


022±0 


019 


0.357±0 


023±0 


020 





482- 


-0 


526 





503 





391±0 


025±0 


023 


0.273±0 


019±0 


018 





526- 


-0 


570 





547 





374±0 


032±0 


028 


0.180±0 


016±0 


014 





570- 


-0 


658 





609 





420±0 


037±0 


036 


0.122±0 


011±0 


011 





658- 


-0 


768 





704 





392±0 


070±0 


049 


0.048±0 


009±0 


006 



Total Observed/Evt. 



1.800±0.016±0.124 



Table 3: Charged kaon fraction and differential cross section per hadronic Z° decay. 
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Xp Range 


<x p > 


/p 


1/N dn p /dxp 


A 


016- 


-0 


022 


A 




A 1 A 

019 


A 




non i n 

029±0 


nn r i n 

005±0 


n i o 

013 


n r r i 1 

9.55±1 


55 


1 A 

±4 


o o 

33 


A 




022- 


A 

-0 


027 


A 




a o r 

025 


A 




n a i in 

041±0 


nno i n 

003±0 


nno 

008 


i n Tn i n 

10.79±0 


O /I 

84 


1 O 

±2 


nn 

09 


A 




027- 


a 

-0 


noo 

033 


A 




030 


A 




A /? | A 

Uo4±U 


nno_i_n 

002±0 


nnr 

005 


1 o et n in 

13.56±0 


A T 

47 


I A 

±0 


no 

98 


A 






-0 


AOO 

038 


A 




no/ 1 

036 


A 




a/ 1 r i n 

065±0 


nno i n 

002±0 


nn a 

004 


11 r /i i n 

11.54±0 


o r 

35 


i n 

±0 


63 


A 




non 

038- 


A 

-0 


C\ A A 

044 


A 




n a 1 

041 


A 




nr* 1 i n 

061±0 


nno i n 

002±0 


nno 

002 


a no i n 

9.03±0 


on 

30 


i n 

±0 


O ET 

25 


A 




C\ A A 

044- 


A 

-0 


049 


A 




n a t 

047 


A 




n/ 1 t i a 

067±0 


nno i n 

002±0 


nno 

002 


o ro i n 

8.52±0 


on 

29 


i n 

±0 


OO 

23 


A 


049- 


A 

-0 


nrr 

055 


A 




n r o 

052 


A 




A /" O 1 A 

062±0 


nno i n 

002±0 


nno 

002 


z 1 oo i n 

6.83±0 


o z 1 

26 


i n 

±0 


OO 

22 


n 

(j 


nrr 

055- 


A 




060 


A 




n c o 

058 


A 




n^o i a 

072±0 


nno i n 

003±0 


nnr 

005 


o r i a 

6.85±0 


OO 

28 


I A 

±0 


A O 

48 


n 
U 


UOU- 


n 
-u 


Uoo 


n 
U 


UOO 


n 
U 


U / 4±U 


nnQ-un 
UUo±U 


UUD 


o. /U±U 


zo 


1 n 
±U 


/IO 
4Z 


n 



Ooo- 


n 

-0 


0/1 


n 



n^n 
009 


n 



n^r i n 

/5±0 


nn /i i n 
004±0 


nn £ 
005 


5.b9±0 


27 


±0 


a n 
40 





U/ 1- 


n 

-0 


Oil 


n 



0/4 


n 



n^7r i n 

/5±0 


nn /i i n 
004±0 


00b 


K no i n 

5.03±0 


27 


±0 


oo 
38 


n 



077- 


n 

-0 


082 





n^7n 

079 





072±0 


nn /I i n 

004±0 


006 


A oo i n 

4.33±0 


27 


±0 


oo 

38 


n 
(J 


082- 


n 
"(J 


Uoo 


A 
(J 


no c 

Uo5 


A 
(J 


no c i n 


nn et i n 

(J(J5±(J 


nn*7 

007 


4.o5±U 


29 


±0 


on 

39 


U 


088- 


-U 


nnn 

uyy 


A 

u 


nno 


n 
U 


nTT_i_n 
Uf f ±U 


nn /I i n 
UU4±U 


nnn 
009 


O P. A i n 


20 


±0 


A 1 

41 


n 
U 


099- 


-U 


i in 

11U 


A 

u 


i n /i 
1U4 


n 
U 


noT_i_n 


nn/^_i_n 
UL)0±U 


mo 
012 


o /i o i n 


23 


±0 


4o 


U 


1 1 n 
110- 


-U 


1 oi 

Izl 


A 

u 


i i r 

llo 


n 
U 


no A i n 


nnT_i_n 


ni k 
015 


O on i n 


25 


±0 


a n 
49 


a 

o 


1 01 

121- 


A 

-0 


1 A O 

143 


A 




1 O 1 

131 


A 




no et i n 

085±0 


nno_i_n 

008±0 


noi 
021 


o oo i n 


21 


±0 


54 


n 

o 


1 /l o 

143- 


n 



a a a 
lo4 


n 



1 £ o 

153 


n 



123±0 


m C i n 

01o±0 


non 
039 


O /I o i n 

z.4z±U 


32 


±0 


/ / 


A 


230- 


A 

-0 


252 


A 




O 1 1 

241 


A 




1 A / J 1 A 

106±0 


nn ^7 i n 

007±0 


ATA 

010 


n *i 0*1 i a 

0.767±0 


048±0 


nT /i 

074 


A 




252- 


A 

-0 


274 


A 




263 


A 




1 1 a i n 

114±0 


nnT i n 

007±0 


ATA 

010 


a nno i n 

0.668±0 


043±0 


n et n 

059 


A 


274- 


A 

-0 


on/ 1 

296 


A 




on c 

285 


A 




-|Ar 1 a 

105±0 


nno i n 

008±0 


nnn 

009 


0.497±0 


036±0 


n a a 

044 


A 


OA/ 1 

296- 


A 

-0 


O 1 O 

318 


A 




O A ^7 

307 


A 




1 A A 1 A 

109±0 


nno i n 

008±0 


nnn 

009 


n a i o in 

0.413±0 


032±0 


no r 

035 





318- 


-0 


351 





334 





099±0 


007±0 


009 


0.296±0 


022±0 


026 





351- 


-0 


384 





366 





098±0 


008±0 


008 


0.219±0 


018±0 


019 





384- 


-0 


417 





400 





105±0 


009±0 


007 


0.175±0 


015±0 


013 





417- 


-0 


450 





432 





104±0 


010±0 


007 


0.134±0 


013±0 


009 





450- 


-0 


482 





465 





103±0 


011±0 


006 


0.101±0 


011±0 


006 





482- 


-0 


526 





503 





095±0 


011±0 


006 


0.066±0 


00^ 


^±0 


004 





526- 


-0 


570 





547 





110±0 


013±0 


006 


0.053±0 


006±0 


003 





570- 


-0 


658 





609 





066±0 


010±0 


006 


0.019±0 


003±0 


002 





658- 


-0 


768 





704 





107±0 


016±0 


007 


0.013±0 


002±0 


001 





768- 


-0 


987 





836 





087±0 


027±0 


012 


0.002±0 


001±0 


000 



Total Observed/Evt. 



0.864±0.015±0.106 



Table 4: Proton fraction and differential cross section per hadronic Z° decay. 
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Xp Range 


<x p > 


1/N dn chg /dxp 





008- 


-0 


010 





009 


509.6±1 


6 


±8 


9 





010- 


-0 


012 





Oil 


481.9±1 


6 


±8 


4 





012- 


-0 


014 





013 


440.9±1 


5 


±7 


7 





014- 


-0 


016 





015 


398.0±1 


4 


±6 


9 





016- 


-0 


022 





019 


334.6±0 


9 


±5 


8 





022- 


-0 


027 





025 


265.2±0 


8 


±4 


6 





027- 


-0 


033 





030 


215.2±0 


7 


±3 


7 





033- 


-0 


038 





036 


178.6±0 


6 


±3 


1 





038- 


-0 


044 





041 


150.0±0 


6 


±2 


6 





044- 


-0 


049 





047 


129.2±0 


5 


±2 


2 





049- 


-0 


055 





052 


111.7±0 


5 


±1 


9 





055- 


-0 


060 





058 


97.2±0 


5 


±1 


7 





060- 


-0 


066 





063 


86.3±0 


4 


±1 


5 





066- 


-0 


071 





069 


77.2±0 


4 


±1 


3 





071- 


-0 


077 





074 


68.7±0 


4 


±1 


2 





077- 


-0 


082 





079 


61.6±0 


4 


±1 








082- 


-0 


088 





085 


56.35±0 


35 


±0 


96 





088- 


-0 


099 





093 


48.53±0 


23 


±0 


83 





099- 


-0 


110 





104 


40.40±0 


21 


±0 


69 





110 


-0 


121 





115 


34.32±0 


20 


±0 


59 





121- 


-0 


143 





131 


27.12±0 


12 


±0 


47 





143- 


-0 


164 





153 


20.35±0 


11 


±0 


35 





164- 


-0 


186 





175 


15.65±0 


09 


±0 


28 





186- 


-0 


208 





197 


12.05±0 


08 


±0 


22 





208- 


-0 


230 





219 


9.50±0 


07 


±0 


17 





230- 


-0 


252 





241 


7.54±0 


07 


±0 


14 





252- 


-0 


274 





263 


6.11±0 


06 


±0 


12 





274- 


-0 


296 





285 


4.969±0 


053±0 


098 





296- 


-0 


318 





307 


3.978±0 


048±0 


081 





318- 


-0 


351 





334 


3.163±0 


035±0 


067 





351- 


-0 


384 





366 


2.367±0 


030±0 


052 





384- 


-0 


417 





400 


1.767±0 


026±0 


041 





417- 


-0 


450 





432 


1.359±0 


023±0 


033 





450- 


-0 


482 





465 


1.028±0 


019±0 


026 





482- 


-0 


526 





503 


0.735±0 


014±0 


020 





526- 


-0 


570 





547 


0.503±0 


012±0 


015 





570- 


-0 


658 





609 


0.300±0 


006±0 


009 





658- 


-0 


768 





704 


0.123±0 


003±0 


004 





768- 


-0 


987 





836 


0.027±0 


001±0 


001 



Table 5: Differential cross section (l/~N)dn c h g /dx p for inclusive charged particle production 
per hadronic Z° decay. The first error is statistical, the second systematic. 
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4.2 Neutral K°/K° and A°/A° Production 

We reconstructed the charged decay modes — > 7r + n~ and A°(A°) — >p7r~(p7r + ) |26 
collectively referred to as V° decays. In order to ensure good invariant mass resolution 
tracks were required to have a minimum transverse momentum of 150 MeV/c with respect 
to the beam direction, at least 40 hits measured in the CDC, and a polar angle satisfying 
|cos#| < 0.8. 

Pairs of oppositely charged tracks satisfying these requirements were combined to form 
V°s if their separation was less than 15 mm at their point of closest approach in 3 dimensions. 
A x 2 fit of the two tracks to a common vertex was performed, and to reject combinatoric 
background we required: the confidence level of the \ 2 t° be greater than 2%; the vertex to be 
separated from the IP by at least 1 mm, and by at least 5ai, where <Ji is the calculated error 
on the separation length of the V^ ; and vertices reconstructed outside the Vertex Detector 
to have at most one VXD hit assigned to each track. 

The two invariant masses m n7T and m p7T were calculated for each V° with, in the latter 
case, the proton (charged pion) mass assigned to the higher- (lower)-momentum track. In 
the plane perpendicular to the beam, the angle between the vector sum of the momenta of 
the two charged tracks and the line joining the IP to the vertex was required to be less than 
both 60 mrad and k ■ (2 + 20 /p± + 5/p 2 ± ) mrad. Here, p± is the component of the vector 
sum momentum transverse to the beam in units of GeV/c and k=1.75 for A°/A° candidates 
and 2.5 for K® candidates. For A°/A° candidates, a minimum vector-sum momentum of 500 
MeV/c was required. 

Note that it is possible for one V° to be considered a candidate for both the K® and 
A /A hypoth eses. Kinematic regions exist where the two hypotheses cannot be distinguished 
without particle identification. In addition there is background from other processes that 
occur away from the IP, most notably 7-conversions into e + e~ pairs. Depending upon the 
type of analysis, such "kinematic-overlaps" may introduce important biases. In this analysis, 
the kinematic-overlap region was removed only when it distorted the relevant invariant mass 
distribution. For the K® analysis, the A°/A° background causes an asymmetric bump in 
the m vix distribution, which complicated the subsequent fitting procedure. A cut on the 7r + 
helicity angle 8* , defined as the angle between the 7r + momentum vector in the rest frame 
and the flight direction, of |cos#*| < 0.8 was used to remove the A , A and 7-conversion 
contamination. 

For the A°/A° analysis, the shape of the K® background depends strongly on momentum. 
Above a V° momentum of a few GeV/c, the K® — > tt + tt~ background is essentially uniform 
in the peak region of the m p7r distribution and no cuts were made to remove the K® overlap. 
At sufficiently low momentum, the K® background becomes asymmetric under the A°/A° 
peak due to detector acceptance; the softer 7r fails to be reconstructed and thus the is not 
found. Therefore, A°/A° candidates with total momentum below 1.8 GeV/c were required to 
have m-xn more than 3a away from the K® mass, where a is the measured resolution on m n7T , 
parameterized as a n7T (p) = 4.6 — 0.27p + 0.21p 2 — O.Olp 3 MeV/c 2 , and p is the V° momentum 
in GeV/c. In order to remove 7 conversions, the proton helicity angle was required to satisfy 
cos#; > -0.95. 

The and m p7r distributions for the remaining candidates are shown in figs. |5| and 
H respectively. The V° candidates were binned in x p , and the resulting invariant mass 
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Figure 5: Invariant mass distribution for all K° s — > 7r + 7r~ candidates. 
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Figure 6: Invariant mass distribution for all A — > p7r and A — > p7r + candidates. 
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distributions were fitted using a sum of signal and background functions. The function 
used for the signal peak was a Gaussian or a sum of two or three Gaussians of common 
center, depending on x p . A single Gaussian was sufficient to describe the K® data in the 
lowest-Xp bin and the A°/A° data in the three lowest-x p bins. However, the mass resolution 
is momentum-dependent and varies substantially over the width of a typical x p bin; two 
Gaussians were sufficient in most cases, with three being needed for both the K® and A°/A° 
data in the highest-x p bin. The relative fractions and nominal widths of the Gaussians 
in the sum were fixed from the MC simulation. The normalization, common center, and a 
resolution scale-factor were free parameters of the fit. The fitted centers were consistent with 
world average mass values |27[], and the fitted scale factor was typically 1.1. The background 
shape used for the K® fits was a quadratic polynomial; for the A / A fits a more complicated 
function was required due to the proximity of the kinematic edge to the signal peak. The 
function P^gim) = a + b{m — m&) + c(l — e d (( m - m A)-o.038)-j wag f ounc [ ^ Q fog adequate in 
Monte Carlo studies, where a,b,c,d were free parameters. 

The efficiencies for reconstructing true and A°/A° decays were calculated, using the 
simulation, by repeating the full selection and analysis on the simulated sample and dividing 
by the number of generated K® or A°/A°. Several checks were performed to verify the MC 
simulation, and thus the V° reconstruction efficiency In particular, the proper lifetimes of 
the and A were measured, yielding values consistent with the respective world averages. 
The simulated reconstruction efficiencies are shown in fig. [j], and were parametrized as func- 
tions of x p . The reconstruction efficiency is limited by the detector acceptance of ~0.67 and 
the charged decay branching fractions of 0.64 for A°/A° and 0.68 for K®. The efficiency at 
high momentum decreases due to finite detector size and two-track detector resolution, and 
the efficiency at low-momentum is limited by the minimum p± and flight distance require- 
ments. The discontinuity in the A°/A° reconstruction efficiency is due to the imposed K® 
mass cut for \ow-x p candidates. 

The differential cross section 1/N dn/dx p per hadronic Z° decay was then calculated in 
each bin by dividing the integrated area under the fitted mass peak by the efficiency, the 
bin width and the number of observed hadronic events corrected for trigger and selection 
efficiency. As is conventional, the K°/K° cross section was obtained by multiplying the 
measured K® cross section by a factor of 2 to account for the undetected K\ component. 
The resulting differential cross sections, including point-to-point systematic errors, discussed 
below, are shown in fig. [12] and listed in table |6|. 

Several sources of systematic uncertainty were investigated for the K® and A°/A° analysis. 
An important contribution to the overall V° spectrum is the track reconstruction efficiency 
of the detector, which was tuned using the world average measured charged multiplicity 
in hadronic Z° decays. We take the ±1.7% normalization uncertainty discussed above (sec. 
4.1) as the uncertainty on our reconstruction efficiency, which corresponds to a normalization 
error on the K°/K° and A°/A° differential cross sections of 3.4%. This uncertainty is 
independent of momentum and is not shown in any of the figures or included in the errors 
listed in table |5|. The momentum-dependent term discussed above and a conservative 50% 
variation of an ad hoc correction |26|] to the simulated efficiency for ^°s that decayed near 
the outer layers of the VXD were also included as systematic uncertainties due to detector 
modelling. 

Each of the cuts used to select V° candidates was varied independently J26J and the 
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Figure 7: The simulated reconstruction efficiencies as a function of x p for (squares) and 
A°/A° (triangles). The charged decay branching ratios are included in the efficiency. The 
discontinuity in the A°/A° reconstruction efficiency at x p = 0.04 is due to the invariant-mass 
cut to remove the low-momentum K® background. 



analysis repeated. For each bin the rms of this set of measurements was calculated and 
assigned as the systematic uncertainty due to modelling of the acceptance. For both the 
K°/K° and the A°/A° candidates, the signal and background shapes used in the fits were 
varied. Single and multiple independent Gaussians, without common centers or fixed widths, 
were used for the signal. Alternative background shapes included constants and polynomials 
of differing orders. In each case the fits were repeated on both data and simulated invariant 
mass distributions and the rms of the resulting differential cross sections was assigned as a 
systematic uncertainty. The MC statistical error on the calculated reconstruction efficiency 
was also assigned as a systematic error. These errors were added in quadrature to give the 
total systematic error. 



4.3 Neutral K*"/K* Q and Production 

We reconstructed the strange vector mesons <ft an d K*°/K*° in the charged decay modes 
<f) — > K + K~ and K*°/K*° — > K ± n T | 28fl . In order to ensure good invariant mass resolution, 



tracks were required to have at least 40 hits measured in the CDC, a track fit quality of 
X 2 /dof< 7, and a polar angle satisfying |cos#| < 0.8. Pairs of oppositely charged tracks 
satisfying these requirements were combined to form neutral candidates if a % 2 fit of the two 
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Neutral V° Production 


x p Range 


< x p > 


1/N dn K ()/dx p 


< x p > 


1/N dn A o/dx p 


0.009-0.011 


0.010 


18.1±1.7 ±2.4 






0.011-0.014 


0.013 


19.1±1.2 ±1.1 






0.014-0.018 


0.016 


20.44±0.91±0.67 


0.015 


2.99±0.45±1.22 


0.018-0.022 


0.020 


21.74±0.85±0.72 


0.020 


3.90±0.42±0.58 


0.022-0.027 


0.025 


20.51±0.70±0.53 


0.025 


4.10±0.30±0.23 


0.027-0.033 


0.030 


17.73±0.55±0.41 


0.030 


3.54±0.23±0.16 


0.033-0.041 


0.037 


16.20±0.46±0.34 


0.037 


3.34±0.20±0.14 


0.041-0.050 


0.045 


13.48±0.38±0.27 


0.045 


2.86±0.14±0.13 


0.050-0.061 


0.055 


11.40±0.31±0.21 


0.055 


2.39±0.11±0.13 


0.061-0.074 


0.067 


10.09±0.27±0.18 


0.067 


2.20±0.10±0.09 


0.074-0.091 


0.082 


8.12±0.23±0.15 


0.082 


1.63±0.08±0.06 


0.091-0.111 


0.100 


6.41±0.20±0.12 


0.100 


1.31±0.08±0.08 


0.111-0.142 


0.126 


4.95±0.16±0.09 


0.125 


0.98±0.06±0.05 


0.142-0.183 


0.161 


3.66±0.16±0.08 


0.160 


0.68±0.05±0.04 


0.183-0.235 


0.206 


2.53±0.17±0.07 


0.205 


0.51±0.05±0.04 


0.235-0.301 


0.262 


1.52±0.08±0.05 


0.262 


0.30±0.04±0.04 


0.301-0.497 


0.371 


0.60±0.05±0.02 


0.368 


0.15±0.02±0.03 


Total Observed/Evt. 


1.90±0.02±0.07 




0.37±0.01±0.02 



Table 6: Measured differential cross sections of neutral K° / K° -mesons and A°/A°-hyperons 
per hadronic Z° decay. A 3.4% normalization uncertainty is included in the systematic errors 
on the observed totals, but not in those on the cross sections. 



tracks to a common vertex converged. The background from long-lived species was rejected 
by requiring the fitted vertex to be within 10 cm or 9<Ji of the IP in three dimensions, and 
within 4 cm or 6<Ji in the plane transverse to the beam direction. The background from 
7-conversions was rejected by assigning the electron mass to both tracks and requiring m ee 
to be greater than 70 MeV/c 2 . 

To reject the high combinatoric background from 7r + 7r~ pairs we used the CRID to 
identify charged kaon candidate tracks. Only liquid (gas) information was used for tracks 
with p < 2.5 (> 3.5) GeV/c, and liquid and gas information was combined for the remaining 
tracks. For this analysis a track was considered "identifiable" if it extrapolated through an 
active region of the appropriate CRID radiator (s); it was considered identified as a kaon 
if the log-likelihood difference between the kaon and pion hypotheses, Ck — Ar, exceeded 
3. These cuts are considerably looser than those used in section 4.1, in order to maximize 
the acceptance for the neutral vector mesons. Efficiencies for identifying selected tracks as 
kaons by this definition were calibrated using the data in a manner similar to that described 
in section 2.2. The K — > K efficiency was found to have a momentum dependence very 
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similar to the 7r — 7r efficiency shown in the upper left plot of fig. |3|, with about 12% lower 
amplitude. There is no dip in the 5-10 GeV/c region since no cut was made against protons. 
The 7r — > K misidentification rate averages 10% and is roughly independent of momentum; 
the p — > K misidentification rate is substantial, especially in the 3-10 GeV/c region, but 
protons constitute only a small part of the combinatoric background. 

A track pair was accepted as a <p ~ > K + K~ candidate if both tracks were identified as 
kaons. A pair was accepted as a K*° — > K + ir~ candidate if one track was identified as a 
kaon and the other was not. Thus a track pair cannot be both a K*°/K*° and candidate. 

The candidates were binned in x p , and the resulting uikk distributions were fitted in a 
manner similar to that described above for the V° candidates. The signal shape was a sum 



of Gaussians of common center; the center was fixed at the world-average mass value f27 
and the amplitude and a resolution scale factor were free parameters. A typical fitted scale 
factor was 1.08. The background shape was parametrized as a threshold term multiplied by 
a slowly decreasing exponential: 



Pbkg(x) = Nx^e ClX+C2x2+C3x3+C4x4+C5x5 (1) 

where x = rrixK — ^ m K, N is an overall normalization factor, and 7 and C1...5 are free 
parameters. Initial values of the background parameters were determined from fits to the 
ttikk distributions for simulated true combinatorial background and for same-sign track pairs 
in the data. The resulting parameters were consistent with each other and the functions 
described the shape of the distribution for candidates in the data in the region away from 
the signal peak. The measured vtikk distributions for the six x p bins are shown in fig. ^, 
along with the results of the fits. 

The case of the K*°/K*° is considerably more complicated due to the natural width of 
the K*° and the presence of many reflections of resonances decaying into 7r + 7r The 
K*°/K*° signal was parametrized using a relativistic Breit-Wigner with the amplitude free 
and the center and width fixed to world-average values |2^j . The background was divided into 
combinatorial and resonant pieces. The combinatorial piece was described by a polynomial 
parametrization similar to that of the <j) but with seven parameters. Parameter values derived 
from fits to simulated combinatorial background and a same-sign data test sample were found 
not to agree with each other or with the opposite-sign data away from the peak, and a search 
over a space of initial values was required in order to find the best fit. 

Knowledge of the resonant contributions to the background is essential, since the K*° is a 
wide state and non-monotonic background variation within its width can lead to systematic 
errors in the measured cross section. We considered four classes of reflections: 

• p° — > 7r + 7r~, K® — ► 7r + 7r~, and uj°,i],ri' — > Ntt, where one of the charged pions is 
misidentified as a K . These backgrounds are large, even after reduction by a factor 
of about 5 by the particle identification. They are particularly important since the 
combination of p and u decays gives rise to a dip in the total background near the 
center of the signal peak, and there is some uncertainty as to the shape of the p 
resonance in Z° decays (see ref. f2~9|). 



7 conversions where one electron is misidentified as a kaon. These are removed effec- 
tively by the m ee cut against 7 conversions noted above. 
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Figure 8: Distributions of invariant mass vtikk for <P candidates in six momentum bins. 
The points with error bars represent the data. The solid curves represent the results of the 
fits described in the text; the dashed curves represent the fitted background component. 
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Figure 9: Simulated relative contributions of the K*°/K*° signal (line) and of various 
resonant backgrounds (dashed lines) to the ttik-k distribution after all analysis cuts. 
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• (j) — > K + K~, where one track is identified as a kaon but the other is not. This 
background is reduced substantially by the requirement that only one of the tracks in 
the pair is identified as a kaon. 

• A — > pir, where the proton is misidentified as a kaon. These are removed effectively 
by the cut against long-lived species noted above. This and the last two categories 
give rise to a more pronounced shoulder in the background just below the signal peak, 
so their removal is quite useful in obtaining a robust fit. 

The shape of the vtik-k distribution for each reflection was parametrized by a smooth 
function fitted to its simulated m^n distribution, and its total production cross section was 
set to the world average value f27j for Z° decays. Figure shows the simulated relative 



contributions from the main resonant backgrounds along with the simulated signal, which 
was scaled to match our measured total cross section (see below). The set of reflection 
functions was added to the combinatorial function to give the total background function. A 
scale factor for each of the four categories of reflections was included as a free parameter in 
the fit to account for possible mismodelling of the misidentification rates; their fitted values 



were consistent with unity. Figure [H] shows the ttlk-k distribution for each momentum bin, 
along with the results of the fits. 

As for the and A°/A° analysis, the (f) and K*°/K*° reconstruction efficiencies were 
determined using the simulation, and are shown in fig. [11]. Differential cross sections were 
calculated in the same way as for the and A°/A°, and the results are shown in fig. [12 
and listed in table |7|. 
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Figure 10: Distributions of invariant mass m^^ for K*°/K*° candidates in six momentum 
bins. The points represent the data. The solid curves represent the results of the fits 
described in the text; the dotted and dashed curves represent the fitted total background 
and combinatoric background components, respectively. 
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Figure 11: The simulated reconstruction efficiencies as a function of x p for (open di- 
amonds) and K*°/K*° (diamonds). The charged decay branching ratios are included in 
the efficiency. The dip in the <fi efficiency at x p ~ 0.13 reflects the dip in the CRID K-7r 
separation at p ~ 2.5 GeV/c (see fig. 3, upper left). 







Neutral Strange 


Meson Production 




x p Range 


<x p > 


1/N d7iK*o/dx p 


Xp Range 


<x p > 


1/N dn^/dxp 


0.018-0.048 


0.033 


4.69 ±0.56 ±0.33 


0.018-0.057 


0.037 


0.744 ±0.074 ±0.048 


0.048-0.088 


0.068 


3.79 ±0.21 ±0.17 


0.057-0.079 


0.068 


0.411 ±0.055 ±0.033 


0.088-0.149 


0.118 


2.23 ±0.13 ±0.14 


0.079-0.175 


0.127 


0.255 ±0.026 ±0.021 


0.149-0.263 


0.206 


1.012±0.056±0.062 


0.175-0.263 


0.215 


0.167 ±0.018 ±0.020 


0.263-0.483 


0.342 


0.343±0.019±0.019 


0.263-0.483 


0.357 


0.0739±0.0068±0.0085 


0.483-1.000 


0.607 


0.051±0.004±0.004 


0.483-1.000 


0.689 


0.0089±0.0015±0.0011 


Total Observed/Evt. 


0.647±0.022±0.029 






0.0985±0.0046±0.0055 



Table 7: Measured differential cross sections of K*°/K*° and mesons per hadronic Z° 
decay. A 3.4% normalization uncertainty is included in the systematic errors on the observed 
totals, but not in those on the cross sections. 
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Systematic uncertainties for this analysis were grouped into efficiency and fit-related 
categories. The dominant contributions to the efficiency category were the uncertainty in the 
track-finding efficiency (see above) and the uncertainty in kaon identification efficiency, for 
which the statistical error on the calibration from the data was used. The total uncertainties 
on the reconstruction efficiencies were 4-6% for K*° /K*° and 6-11% for 0, depending on 
momentum. 

In the case of the 0, fitting systematics were evaluated by varying the signal shape as 
in the V° analysis. In addition, fits were performed with the signal center shifted by plus 
and minus the error on the world-average mass value. The effect of background fluctuations 
was evaluated by taking the largest variation in the result over a set of fits done with the 
background shape parameters Cj fixed to all combinations of their fitted values ±ler. The 
total fitting uncertainties were 2-8%. 

In the case of the K*°/K*°, we considered the same variations, as well as variation of 
the signal width by ±lcx from the world-average value and several variations of the resonant 
background. Fits were performed with the misidentification scale factors fixed to their fitted 
values ±50% for the tttc category and ±15% for the others, corresponding to roughly twice 
the error on our measured misidentification rates. All 16 combinations were considered, and 
the largest variation taken as a systematic error. The cross section for production of each 
resonance was varied by the error on the world-average value. The sizes of the p and u 
contributions were varied in all four combinations of ±30% and ±10%, respectively, and 



the largest variation was taken as a systematic error. Following [[29[] an error due to the 
uncertainty in the p° lineshape was evaluated by shifting the p reflection function down by 
40 MeV/c 2 . The total fitting uncertainties were 2-6%. 

4.4 Hadron Production in Inclusive Hadronic Z° Decays 

Our measured differential cross sections per hadronic event of the seven hadron species are 



shown as a function of x p in fig. 12, along with that of inclusive charged particles. At low 



x p pions are seen to dominate the hadrons produced in hadronic Z° decays. For example, at 
x p ~ 0.03, pseudoscalar K ± and K°/K° are produced at a rate about ten times lower than 
pions, vector K*° are suppressed by an additional factor of ~4, and the doubly strange vector 
by another factor of ~12. The most commonly produced baryons, protons, are suppressed 
by a factor of ~25 relative to pions, and the strange baryon A°/A° by an additional factor 
of ~3. 

These results are in general consistent with previous measurements from experiments at 
LEP [|7| , provided that the point-to-point correlations in the systematic errors are taken into 
account. However, although our proton differential cross section for x p > 0.35 is consistent 



with that measured by ALEPH | 24fl , it is not consistent with that measured by OPAL [^3 

We compared our results with the predictions of the JETSET 7.4, UCLA 4.1 and HER- 
WIG 5.8 event generators described in section 1, using in all cases the default parameters. 



Figures O and 114] show the charged fractions and the neutral differential cross sections, 



respectively, along with the predictions of these three models. The momentum dependence 
for each of the seven hadron species is reproduced qualitatively by all models. For momenta 
below about 1.5 GeV/c, all models overestimate the kaon fraction significantly and all ex- 
cept UCLA underestimate the pion fraction by about 2a (taking into account the correlation 
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Figure 12: Differential cross sections per hadronic Z° decay per unit x p for inclusive charged 
particles (open circles), ir^ (circles), K ± (squares), K°/K° (open squares), K*°/K*° (dia- 
monds), (open diamonds), p/p (triangles), and A°/A° (open triangles). The baryon and 
all-charged differential cross sections have been scaled by 0.04 and 1.1, respectively, for clar- 
ity. The error bars represent statistical and systematic errors added in quadrature. Each 
point is plotted at the average x p value of reconstructed particles in that bin (see tables 2-7). 
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Figure 13: Comparison of our measured charged hadron fractions (symbols) with the 
predictions of the JETSET (dashed lines), UCLA (solid lines) and HERWIG (dotted lines) 
fragmentation models. 



30 




0.01 0.1 1 

4-98 y 
8400 A9 A p 



Figure 14: Comparison of our measured neutral hadron differential cross sections with the 
predictions of three fragmentation models. 



31 



in the experimental errors). In the 5-10 GeV/c range UCLA and HERWIG overestimate 
the pion fraction by 2— 3a. For p > 10 GeV/c, JETSET overestimates the proton fraction, 
but describes the momentum dependence. In this momentum region, HERWIG and UCLA 
predict a momentum dependence in the proton fraction that is inconsistent with the data. 

In the case of K°/K°, all models describe the data well at high x p , but overestimate 
the cross section at low x p by as much as 50%. A similar excess was seen in the charged 
kaon fraction (see fig |T3]) . In the case of A°/A°, JETSET and UCLA describe the data well 
except for a 10% shortfall near x p = 0.02. HERWIG describes the data well except for the 
lowest and highest x p points, where it overestimates the production. The structure in the 
HERWIG prediction at very high x p is similar to that seen in the proton fraction, and is also 
visible to varying degrees in the predictions for the neutral strange mesons. In the case of 
K*°/K*°, JETSET is high by a roughly constant factor of 1.5 across the x p range; HERWIG 
and UCLA reproduce the data except at the lowest x p point. In the case of 0, JETSET is 
high by a factor of two over all x p , UCLA is high for x p > 0.06, and HERWIG describes the 
data except at the highest x p point. 



5 Flavor-Dependent Analysis 

The analyses described above were repeated on the light-, c- and 6-tagged event samples 
described in section 3, to yield differential cross sections R^ a9 for each hadron species h in 
each tagged sample. True differential cross sections R™ in events of the three flavor types, 
k,m — I, c, b, representing events of the types Z° — ► uu,dd,ss, Z° — * cc, and Z° — ► bb, 
respectively, were extracted by solving for each species h the relations: 

Rh = ^ m e mk F m ■ (2) 

Here, F m is the fraction of hadronic Z° decays of flavor type m, taken from the Standard 
Model, e m k is the event tagging efficiency matrix (see table |l|), and B^ k represents the 
momentum-dependent bias of tag k toward selecting events of flavor m that contain hadrons 
of species h. Ideally all biases would be unity in this formulation. The biases were calculated 
from the MC simulation as B^ k = (n^ ktag /N mthtag )/ (n^/N m ), where N m (n£J is the number 
of simulated events (hadrons of species h in events) of true flavor m and N m ^ tag (n^ ktag ) 
is the number of (/i-hadrons in) those events that are tagged as flavor k. The diagonal bias 
values |^0|, [28| are within a few percent of unity for the charged hadrons, <fi and K*°, 
reflecting a small multiplicity dependence of the flavor tags. They deviate by as much as 
10% from unity for the K°/K° and A /A , since some tracks from V° decays are included 
in the tagging track sample and have large impact parameter. The off-diagonal bias values 
deviate from unity by a larger amount, but these have little effect on the unfolded results. 

The resulting differential cross sections are listed in tables |S|-|l"4]. The systematic errors 
listed are only those relevant for the comparison of different flavors, namely those due to 
uncertainties in the unfolding procedure; the systematic errors given in the preceding section 
are also applicable, but are common to all three flavor categories. The flavor unfolding 
systematic errors were evaluated by varying each element of the event tagging efficiency 
matrix en by ±0.01 [fjQfl , varying the heavy quark production fractions Rb and R c by the errors 
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X 
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Ratios 
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cc 
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Table 8: Measured differential cross sections (l/N)dn 7r ±/da;p for the production of charged 
pions per Z° decay into light (u, d, s), c and b primary flavors. The errors are the sum in 
quadrature of statistical errors and those systematic uncertainties arising from the unfolding 
procedure. Systematic errors common to the three flavors are not included. The < x p > 
values for the three flavor samples are consistent in each bin, and have been averaged. 
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Table 9: Differential cross sections for the production of mesons per Z° decay into light, 
c and b primary flavors. 



Xp 




K*°/K*° Production Cross Sections 


Ratios 


Range 


<x p > 


uu, dd, ss 


cc 


bb 


cuds 


b.uds 


0.018-0.048 
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5.2±1.3 
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0.81±0.12 
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1.75±0.40 
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0.342 


0.345±0.042 


0.29±0.20 
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0.85±0.58 


1.16±0.27 


0.483-1.000 


0.607 


0.076±0.010 


0.026±0.034 


0.012±0.009 


0.36±0.45 


0.15±0.11 



Table 10: Differential cross sections for the production of K*°/K*° mesons per Z° decay 
into light, c and b primary flavors. 
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Table 11: Differential cross sections for the production of p/p per Z° decay into light, c 
and b primary flavors. 
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A°/A° Production Cross Sections 
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Table 12: Differential cross sections for the production of A°/A° per Z° decay into light, c 
and b primary flavors. 
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Table 13: Differential cross sections for the production of K°/K° mesons per Z° decay into 
light, c and b primary flavors. 
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Table 14: Differential cross sections for the production of mesons per Z° decay into light, 
c and b primary flavors. 
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on their respective world averages, and varying each diagonal bias value by the larger 
of ±0.005 and ±20% of its difference from unity. Since the lepton background is strongly 
flavor- dependent, the photon conversion rate in the simulation was varied by ±15%, and the 
simulated rates of lepton production from other sources in light-, c, and 6-flavor events were 
varied by ±50%, ±10% and ±5%, respectively. The unfolding systematic errors are typically 
small compared with the statistical errors, and are dominated by the variation in the bias. 

In figure 15 we show the differential cross sections for the seven hadron species in light- 
flavor Z° decays. Qualitatively these are similar to those in flavor-inclusive decays (fig. |T2"D , 
although all differential cross sections are larger at high x p in light flavor events. The 
same general features of tt-K and p-A° convergence at high x p are visible, and the relative 
suppressions of hadron species with respect to one another are similar in magnitude and 
momentum dependence. 

Also shown in fig. [15| are the predictions of the three simulation programs. All models 
reproduce the shape of each differential cross section qualitatively. The JETSET prediction 
for charged pions is smaller than the data in the range x p < 0.015, and those for the 
pseudoscalar kaons are larger than the data for 0.015 < x p < 0.03; those for the vector 
mesons and protons reproduce the x p dependence but show a larger normalization than 
the data. These differences were all seen in the flavor- inclusive results (figs. [13], [14D, and 
we can now conclude that they all indicate problems with the modelling of light-flavor 
fragmentation, and cannot be due entirely to mismodelling of heavy hadron production and 
decay. The HERWIG prediction for pseudoscalar kaons is also larger than the data at low x p 
and is slightly smaller than the data in the range 0.15 < x p < 0.25. For all hadron species the 
HERWIG prediction is larger than the data for x p > 0.4, showing a characteristic shoulder 
structure. The UCLA predictions for the baryons and the vector mesons show a similar but 
less pronounced structure that is inconsistent with the proton and K*°/K*° data. Otherwise 
UCLA reproduces the data except for pseudoscalar kaons in the range 0.015 < x p < 0.03. 

In fig. [16| and tables [8|-[14] we give the ratios of production in 6-flavor to light-flavor events 
for the seven species. The systematic errors on the hadron reconstruction and identification 
largely cancel in these ratios, and the total errors are predominantly statistical. There is 
higher production of charged pions in 6-flavor events than in light-flavor events at low x p , 
with the ratio rising with x p for 0.008 < x p < 0.03 to a plateau value of about 1.25. The 
production of both charged and neutral kaons is approximately equal in the two samples for 
x p < 0.03, but the relative production in 6-flavor events then increases with x p , peaking at a 
value of about 1.7 at x p ~ 0.09. The errors on the K*°/K*° and ratios are large, but the 
data are consistent with behavior similar to that of the pseudoscalar kaon ratios. There is 
approximately equal production of baryons in 6-flavor and light-flavor events for x p < 0.15. 
The production of pions and pseudoscalar kaons in 6-flavor events falls rapidly with x p for 
x p > 0.1 relative to that in light-flavor events. The relative production of the vector mesons 
and protons also falls at high x p . 

These features are consistent with expectations based on the known properties of e + e~ — > 
bb events, namely that a large fraction of the event energy (on average about 70% |7[) is 
carried by the leading B- and B-hadrons, leaving little energy available to produce high 
momentum fragmentation hadrons. The B hadrons decay into a large number of lighter 
particles, including on average 5.5 stable charged hadrons P7| , which are expected to pop- 
ulate primarily the region 0.02 < x p < 0.2. Also shown in fig. [16] are the predictions of 
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Figure 15: Identified hadron differential cross sections in light-flavor events. Also shown 
are the predictions of the three fragmentation models; the prediction of each model for K ± 
is similar to that for K°/K°, and the two have been averaged. 
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Figure 16: Ratios of production of each hadron species in 6-flavor events to that in light- 
flavor events (left) and in c-flavor:light-flavor events (right). Also shown are the predictions 
of the three fragmentation models; for each model, the predictions for K ± and K°/K° were 
averaged, as were those for p/p and A°/A°. The model predictions for are not shown, 
but have the same x v dependence as the corresponding prediction for K*°/K*°, with a peak 
value typically higher by 40%. 
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the three fragmentation models, all of which reproduce these features qualitatively, although 
HERWIG overestimates the ratio for pions in the range x p < 0.05 and that for kaons for 
x p < 0.3. The values of these ratios depend on details of the B and D hadron energy spectra 



and decay properties, and so provide information complementary to that in fig. [15]. How- 
ever, in drawing conclusions regarding heavy flavor modelling from these ratios, one must 
consider how well the model in question reproduces the light flavor results. For example, the 
HERWIG prediction for pion (kaon) production in light-flavor events (fig. 15) is consistent 
with (higher than) the data for x p < 0.05, so it is safe to conclude from fig. [16] that HERWIG 
mismodels pion and kaon production from B decays in this region. However the fact that 
the HERWIG ratio for kaons is high in the region 0.1 < x p < 0.3 is due at least in part to 
the low HERWIG prediction for kaon production in light-flavor events in that region. 

In fig. we also show the ratios of production in c-flavor to light-flavor events for the 
seven species. The errors are larger than for the b:uds comparison and x p bins have been 
combined in some cases for clarity. Similar qualitative features are observed: there is higher 
kaon production in c-flavor events than in light-flavor events at x p ~ 0.1; pion production is 
slightly higher in c-flavor than in light-flavor events for x p < 0.03, then decreases slowly with 
x p ; both pion and kaon production appear to fall rapidly with x p for x p > 0.3, a somewhat 
higher value than the corresponding kuds ratios. These features are expected since c-jets 
produce a charmed hadron with on average about half 0] the beam energy, a lower fraction 
than 5-hadrons, which leaves more energy available for fragmentation hadrons than in fe-jets. 
The charmed hadron decay products often include a kaon carrying a large fraction of the 
charmed hadron momentum, and there are fewer additional charged pions than in B hadron 



decays. Also shown in fig. [16] are the cuds ratios predicted by the three fragmentation 



models. All models are consistent with the data, except that HERWIG overestimates the 
pion ratio for 0.03 < x p < 0.15. 



6 Comparison with QCD Predictions 

We tested the predictions of Gribov and Lipatov, that, in the limit x p — > 1, the momentum 
distribution for primary leading hadrons be (1 — x p ) n , with n = 2 for mesons and n = 3 
for baryons. Since this test benefits from more bins at high x p , we considered only the 
charged hadrons. The cross sections measured in light flavor events provide in principle a 
better test than those measured in flavor-inclusive events, since c- and 6-flavor events cannot 
contain primary leading pions, kaons or protons. However, we have just shown that the 
contributions from c- and 6-flavor events are small for x p greater than about 0.5; since we 
have better statistics for flavor-inclusive events we performed the test on this data set, as 
well as on the light-flavor data. We are limited to x p < 0.77 for the charged pions and 
kaons, but for the flavor-inclusive analysis of protons we have an additional bin, obtained 
from a 2- hypothesis analysis (see sec. 4.1) that also yielded the sum of meson cross sections 
(tt +K ). We also considered this meson sum at all momenta, which has smaller statistical 
errors than the sum of the individual and K cross sections. 

Figure [17] shows the 7r ± , K^, p and (7r ± + K^) differential cross sections as functions of 



(1 — x p ) in flavor-inclusive Z° decays. Fits of the function f(x) = A(l — x p ) n , with the value 
of n fixed to 2 (3 for protons), were performed to the first m data points and the resulting 
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Figure 17: Measured differential cross sections in flavor-inclusive Z° decays as a function 
of (1 — x p ), along with the results of polynomial fits, described in the text, to the data in 
the 2, 4 and 6 leftmost bins. Each fitted polynomial has been integrated over each bin and 
is shown as a histogram. 
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Figure 18: Distribution of £ = ln(l/x p ) for charged kaons in flavor-inclusive Z° decays. 
The solid and dashed lines indicated the results of fits of the Gaussian and distorted Gaus- 
sian approximations of MLLA QCD described in the text. The dotted lines indicate the 
continuations of the fitted Gaussian function. 

fitted distributions for m = 2, 4, 6 are shown in the figure. In all cases the fit quality is good 
for m = 2, but worsens with increasing m. The maximum number of bins for which the 
confidence level of the \ 2 °f the fit exceeded 0.01 was 3 for 7r ± and K ± , 6 for p/p, and 2 for 
the meson sum (7r ± + K^). 

Using this criterion, the theoretical prediction is consistent with our combined meson 
data for (1 — x p ) < 0.34, with our pion and kaon data for (1 — x p ) < 0.47, and with our 
proton data for (1 — x p ) < 0.57. A similar analysis of the light-flavor sample (not shown) 
yielded similar results; the prediction is consistent with our pion, kaon and combined meson 
data for (1 — x p ) < 0.53, and with our proton data for (1 — x p ) < 0.62. 

In order to test the predictions of QCD in the Modified Leading Logarithm Approxi- 
mation (MLLA) combined with the ansatz of Local Parton-Hadron Duality (LPHD), we 
converted our measurements into differential cross sections in the variable £ = ln(l/x p ). 



Figure ^ shows our measured differential cross section as a function of £ for the charged 
kaons. Also shown are the results of fits to a simple Gaussian, and a distorted Gaussian 
including skewness and kurtosis terms. The Gaussian fit was performed over a £ range of 
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-- A4 £ = 3n(1/x p ) ^ = 8n(1/xp) 

Figure 19: Distributions of £ for the seven hadron species in flavor- inclusive hadronic Z° 
decays (points), along with the results of Gaussian fits (solid lines) to the data over a range 
of approximately ±1 unit about the peak. 
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Figure 20: Peak positions from fits to the £ distributions in flavor-inclusive and light- 
flavor hadronic Z° decays. Also shown are averages of similar flavor-inclusive results from 
experiments at LEP. The line is the result of an ad hoc exponential fit to our light-flavor 
data. 



width 2 units positioned near the maximum of the distribution. The fitted peak position 
£* was found to be independent of the exact position of this range within statistical errors, 
and the solid line in fig. |18| represents the result of a fit over a range centered on this peak 
position. A good fit quality was obtained; the two points above this £ range could be added 
to the fit, as could the first two points below the range, before the y 2 began to increase 
rapidly, indicating that the Gaussian approximation is consistent with our data over a range 
of approximately ±f.3 units of £ around the peak position. The distorted Gaussian function 
is able to describe the data over the full measured range of £, as indicated by the dashed line 
in fig. [H| however the distortion terms grow rapidly as points outside the range described 
by the simple Gaussian are added. 

Similar results were obtained for the other hadron species. Their ^-distributions are 
shown in fig. We fitted a simple Gaussian over a £ range of approximately ±1 unit 



centered on the maximum of each distribution in order to measure the peak position £* for 
each hadron species. Systematic errors on this measurement were evaluated by varying the fit 
range and by refitting with each source of correlated experimental systematic error considered 
coherently in turn. Good fit qualities were obtained when the correlated systematic errors 
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Table 15: Peak positions £* from Gaussian fits to the £ distributions for each hadron species 
measured in flavor-inclusive and flavor-specific hadronic Z° decays. The errors are the sum 
in quadrature of statistical and systematic components. 



were taken into account. The peak positions are given in table 15 and shown as a function 
of hadron mass in fig. [2(| along with averages of similar measurements from experiments 
at LEP [0, with which they are consistent. The distribution for pions peaks at a higher 
£ value than the those of the other hadron species, but otherwise there is no monotonic 
mass-dependence. 

As discussed in section 1, the MLLA QCD+LPHD prediction is valid for primary frag- 
mentation particles, whereas experiments so far have measured samples that include decay 
products of an unknown mix of resonances as well as of heavy hadrons. This mix may affect 
measured £* values differently for different hadron species. It is of interest to try to resolve 
this question experimentally, and we have therefore applied the same analysis to the three 
primary event flavor categories discussed in the previous section. We expect the light flavor 
events to be less affected by decay products, as D- and £>-hadron decays are excluded. 

The Gaussian function provides an acceptable description of the £ distribution for all 
hadron species in events of each flavor within about ±1 unit of the peak (not shown), and 
the fitted peak positions are listed in table [15|. For the K*°/K*° and in c-flavor events, 
the limited sample size did not allow a reasonable systematic error evaluation, so they are 
omitted. 

The £* values measured in 6-flavor events are significantly different from those measured 
in light-flavor events for and K*°/K*°; the difference is 1.5a for and K°/K°. For the 
other hadron species the £* values measured in events of all three flavors are consistent. The 
£* values measured in light-flavor events differ significantly from those measured in flavor- 
inclusive events for and p/p. The light-flavor values are also shown in fig. |2(| The 



result of an ad hoc exponential fit to the light-flavor data is shown in fig. |20] as a reference 
trajectory, and the light-flavor data are seen to lie closer to a monotonic trajectory than the 
flavor-inclusive data. 
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Total Cross Sections per Event of Flavor 
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Table f6: Corrected total cross sections per hadronic Z° decay, and per decay into light, 
c or b primary flavor. Differences between the total cross sections for c- and light-flavor 
and b- and light-flavor events. All errors are the sum in quadrature of experimental and 
extrapolation uncertainties 



7 Total Production Cross Sections 

We have integrated our differential cross sections over their respective measurement ranges, 
taking into account the bin-to-bin correlations in the systematic errors. These integrated 
cross sections per event are listed in tables |2|-|7|; the errors are dominated by overall normal- 
ization uncertainties corresponding to the uncertainty in our track reconstruction efficiency. 
In order to quote total cross sections, we must extrapolate into the unmeasured regions of x p , 
and we have done this using the three MC models discussed above. From the hadrons of each 
species generated using each of these models, we calculated the fraction that were generated 
with x p in the range of our measurement. For each hadron species the three fractions were 
found to be similar, with the UCLA (HERWIG) fraction being typically 1% larger (1-2% 
smaller) than the JETSET fraction. The average of the three accepted fractions ranged from 
0.812 for K to 0.945 for K°/K°. Each integrated measured cross section was divided by 
the corresponding average fraction, and an uncertainty of ±0.01 (±0.015) was assigned to 
the average fraction for ir ± , K ± , K°/K°, p/p and A°/A° (K*°/K*° and 0), corresponding 
to a typical rms among the three predictions. The corrected total cross sections are shown 
in table [16], and were found to be consistent with an average of similar measurements from 
experiments at LEP 0. 

As a cross check, we fitted the distorted Gaussian function described in section 6 to the 
£ distribution for each hadron species, and calculated the fraction of the area under the 
fitted curve that was within the range of our measurement. An uncertainty was assigned 
corresponding to the largest variation obtained by varying the fitted parameter values by 
all combinations of ±lcr and —la. The resulting fractions are consistent with those ob- 
tained using the fragmentation models, giving confidence in both the central values and the 
uncertainties assigned. 

We applied the same procedure to our measurements for the three flavor categories. 
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The three simulations were found to give similar flavor dependences, with the accepted 
fraction in b (c) events typically 0.02 (0.01) larger than that in light-flavor events. The 
resulting total cross sections are listed in table ITB along with differences between flavors, 
for which some of the systematic errors cancel. We observe roughly 15% more pseudoscalar 
mesons in 6-flavor events than in light-flavor events, and the respective sums of the charged 
hadron differences are consistent with our previous measurement [3lJ of the differences in 
total charged multiplicity between light-, c- and 6-flavor events. All other differences are 
consistent with zero. 



8 Leading Particle Effects 



We extended these studies to look for differences between particle and antiparticle production 
in light quark (as opposed to antiquark) jets, in order to address the question of whether 
e.g. a primary u-initiated jet contains more hadrons that contain a valence w-quark (e.g. 
7r + , K + , p, A ) than hadrons that do not (e.g. n~ , K~, p, A ). To this end we used the 
light quark- and antiquark-tagged hemispheres described in section 3. 
We measured the differential cross sections per light quark jet 
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where: q and q represent light-flavor quark and antiquark jets respectively; N evts is the total 
number of events in the sample; h represents any of the identified hadron species 7r~, K~, 
K , p, or A , and h indicates the corresponding antihadron. Then, for example, N(q — > h) 
is the number of hadrons of species h in light quark jets. This formulation assumes CP 
symmetry, i.e. N(q — ► h) — N(q — > h), which was found to be satisfied in the data in all 
cases. 

The charged hadron fractions analysis was repeated on the sample of positively charged 
tracks in the quark-tagged jets and negatively charged tracks in the antiquark-tagged jets, 
yielding measured values of R^+, R q K +, and R^ in the tagged samples. The same procedure 
applied to the remaining tracks yielded , R q K - , and Rp. The K* / K*° and A / A analyses 



R K *o, R\ and 



were applied similarly to the quark- and antiquark-tagged jets to yield R—, 
R\. 

The light-tagged event sample contains a residual heavy flavor background of 12% cc and 
3% bb events. The decays of the leading heavy hadrons in simulated heavy flavor background 
events give rise to substantial differences between hadron and antihadron production in the 
quark-tagged sample over the entire x p range. It is essential to understand this contribution, 
which is typically 15% of the observed hadrons for x p < 0.5 and decreases at higher x p 



(see fig. [16]) • The simulated contribution to each cross section was applied as a correction, 
yielding differential cross sections per light-quark-tagged jet. 

For each hadron species, differential cross sections in light quark jets were then extracted 
by correcting for the light-tag bias (see sec. 5) and unfolding for the effective quark (vs. 
antiquark) purity. The purity was estimated from the simulation to be 0.76 for the A / A 
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X 


V 




7r + and it 


Production in u,d,s Jets 


Range 


<x p > 


7T + 


7T 


n 


0.016- 


-0.022 


0.019 


140.9±2.5 


139.0±2.6 


-0.007±0.016 


0.022- 


-0.033 


0.027 


98.2±1.5 


96.7±1.4 


-0.007±0.014 


0.033- 


-0.044 


0.038 


62.8±1.3 


63.6±1.3 


0.007±0.019 


0.044- 


-0.055 


0.049 


44.2±1.4 


A A /~\ 1 -I A 

44.9±1.4 


0.007±0.029 


0.055- 


-0.066 


0.060 


33.4±1.1 


33.2±1.1 


-0.003±0.030 


0.066- 


-0.077 


0.071 


25.79±0.82 


27.16±0.82 


0.026±0.028 


0.077- 


-0.088 


0.082 


21.66±0.71 


22.34±0.71 


0.016±0.029 


0.088- 


-0.099 


0.093 


17.17±0.62 


18.40±0.63 


0.034±0.032 


0.099- 


-0.110 


0.104 


14.45±0.57 


14.52±0.57 


0.003±0.036 


0.110- 


-0.121 


0.115 


11.44±0.50 


12.84±0.52 


0.057±0.038 


0.121- 


-0.143 


0.131 


9.32±0.32 


9.61±0.32 


0.015±0.031 


0.143- 


-0.164 


0.153 


7.21±0.28 


7.39±0.28 


0.012±0.035 


0.164- 


-0.186 


0.175 


5.40±0.24 


5.49±0.25 


0.008±0.041 


0.186- 


-0.208 


0.197 


4.30±0.21 


4.44±0.22 


0.016±0.045 


0.208- 


-0.230 


0.219 


3.14±0.19 


3 30±0 19 


026±0 053 


0.230- 


-0.274 


0.251 


2.37±0.12 


2.59±0.12 


0.043±0.043 


0.274- 


-0.318 


0.295 


1.398±0.091 


1.687±0.097 


0.093±0.055 


0.318- 


-0.384 


0.348 


0.972±0.061 


0.996±0.064 


0.012±0.057 


0.384- 


-0.471 


0.423 


0.456±0.040 


0.504±0.042 


0.050±0.077 


0.471- 


-0.603 


0.527 


0.180±0.025 


0.210±0.026 


0.08±0.12 


0.603- 


-0.768 


0.668 


0.065±0.019 


0.089±0.021 


0.16±0.23 



Table 17: Differential cross sections for the production of positive and negative pions in 
light (u, d and s) quark jets from hadronic Z° decays, along with the normalized difference 
D n - between the two. The errors are the sum in quadrature of statistical errors and those 
systematic errors arising from the light quark tagging and unfolding procedure. 



and 0.72 for the charged hadrons and K*°/K*°, the latter value reflecting the cutoff in 
acceptance of the CRID at | cos^l = 0.68. 

The measured differential cross sections per light quark jet are listed in tables |I~7]-|2l"| for 
the five measured hadron species that are not self-conjugate. As for the flavor dependent 
results (sec. 5), the error given is the sum in quadrature of the statistical error and those 
systematic errors arising from the tagging and correction procedures. The latter include 
variation of the event tagging efficiencies and biases as described in section 5, variation of 
the electroweak parameters -R&, R c , Af, and A c by the errors on their respective world average 
values |27|, and variation of the effective quark purity by ±0.015 to cover the uncertainty in 



the electron beam polarization and statistical error on the simulated purity. The systematic 
errors are small compared with the statistical errors, and are typically dominated by the 
uncertainty on the effective quark purity. These results supersede those in our previous 
publication M. 
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K*° and K*° Production in u,d,s Jets 


Range 


< x p > 


K*° 


K*° 


D R . 


0.018-0.048 


0.033 


2.50±0.94 


2.69±0.95 


0.04±0.29 


0.048-0.088 


0.068 


1.64±0.36 


2.40±0.38 


0.18±0.14 


0.088-0.149 


0.118 


1.11±0.22 


0.88±0.22 


-0.11±0.17 


0.149-0.263 


0.206 


0.318±0.087 


0.447±0.095 


0.17±0.19 


0.263-0.483 


0.342 


0.053±0.033 


0.264±0.042 


0.67±0.18 


0.483-1.000 


0.607 


0.022±0.012 


0.100±0.015 


0.64±0.16 



Table 18: Differential cross sections for the production of K and K mesons in light 
quark jets, along with their normalized difference. 





Xp 

Rang 


e 


< x p > 


K + and K~ 
K+ 


Production in 
K~ 


u,d 


,s Jets 
D K - 





016- 


-0 


022 





019 


8.3±1 


1 


14.8±1 


3 





28±0 


09 





022- 


-0 


033 





027 


9.27±0 


69 


8.14±0 


68 


-0 


06±0 


07 





033- 


-0 


044 





038 


8.05±0 


68 


7.70±0 


68 


-0 


02±0 


08 





044- 


-0 


055 





049 


8.03±0 


81 


7.59±0 


81 


-0 


03±0 


09 





055- 


-0 


066 





060 


3.75±0 


74 


6.27±0 


79 





25±0 


14 





066- 


-0 


088 





077 


3.44±0 


45 


3.90±0 


47 





06±0 


11 





088- 


-0 


121 





101 


3.09±0 


41 


2.73±0 


42 


-0 


06±0 


13 





208- 


-0 


230 





219 


0.99±0 


18 


1.36±0 


19 





15±0 


14 





230- 


-0 


274 





251 


0.595±0 


091 


1.120±0 


099 





31±0 


10 





274- 


-0 


318 





295 


0.383±0 


072 


0.895±0 


081 





40±0 


11 





318- 


-0 


384 





348 


0.260±0 


049 


0.665±0 


055 





44±0 


10 





384- 


-0 


471 





423 


0.163±0 


034 


0.427±0 


039 





45±0 


11 





471- 


-0 


603 





527 


0.091±0 


023 


0.219±0 


026 





42±0 


14 





603- 


-0 


768 





668 


-0.007±0 


017 


0.120±0 


022 


1 


12±0 


28 



Table 19: Differential cross sections for the production of positive and negative kaons in 
light quark jets from hadronic Z° decays, along with their normalized difference. 
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X 


V 








p and p Production in u,d,s 


Jets 






Rang 


e 


< x p > 


P 




P 






D p 









-u 


033 


U 




7.1±1 


1 


4.7±1 


4 





20±0 


21 





033- 


-0 


044 





038 


5.76±0 


52 


4.83±0 


51 





09±0 


09 





044- 


-0 


055 





049 


4.10±0 


44 


4.07±0 


44 





00±0 


10 





055- 


-0 


066 





060 


3.65±0 


44 


3.20±0 


44 





07±0 


12 


U 


066- 


-0 


Uoo 





077 


2.69±0 


30 


2.31±0 


30 





no i n 
Uo±U 


1 1 
1 1 





088- 


-0 


121 





101 


1.82±0 


29 


1.99±0 


30 


-0 


04±0 


14 





230- 


-0 


274 





251 


0.618±0 


078 


0.292±0 


072 





36±0 


15 





274- 


-0 


318 





295 


0.387±0 


056 


0.157±0 


053 





42±0 


18 





318- 


-0 


384 





348 


0.257±0 


035 


0.099±0 


033 





44±0 


18 





384- 


-0 


471 





423 


0.117±0 


020 


0.076±0 


019 





21±0 


19 





471- 


-0 


603 





527 


0.070±0 


010 


0.025±0 


009 





47±0 


19 





603- 


-0 


768 





668 


0.018±0 


004 


0.001±0 


004 





85±0 


42 



Table 20: Differential cross sections for the production of protons and antiprotons in light 
quark jets, along with their normalized difference. 



Xp 




A and A Production in 


u,d,s Jets 


Range 


< x p > 


A 


A 


D A » 


0.010-0.030 


0.022 


0.65±0.16 


1.05±0.17 


-0.23±0.18 


0.030-0.050 


0.040 


0.86±0.13 


0.91±0.13 


-0.03±0.14 


0.050-0.070 


0.060 


0.529±0.084 


0.555±0.084 


-0.02±0.14 


0.070-0.100 


0.083 


0.303±0.057 


0.468±0.060 


-0.21±0.14 


0.100-0.140 


0.118 


0.301±0.053 


0.319±0.054 


-0.03±0.16 


0.140-0.180 


0.158 


0.190±0.048 


0.157±0.047 


0.09±0.25 


0.180-0.300 


0.227 


0.171±0.034 


0.098±0.032 


0.27±0.23 


0.300-0.500 


0.368 


0.090±0.022 


0.013±0.019 


0.75±0.37 



Table 21: Differential cross sections for the production of A and A hyperons in light quark 
jets, along with their normalized difference. 
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It is convenient to show these results in the form of the difference between hadron h and 
antihadron h production normalized by the sum: 

D h = jjq - R 9 • ( 5 ) 

The common systematic errors cancel explicitly in this variable, which is shown for each 



hadron species in fig. A value of zero corresponds to equal production of hadron and 
antihadron, whereas a value of +(-)l corresponds to complete dominance of (anti)particle 
production. In each case the difference is consistent with zero at low x p . For charged pions it 
is also consistent with zero at high x p , but for the other hadrons there are significant positive 
differences that appear to increase with increasing x p . 

The results for the baryons (fig. pl|a,b) afford the most straightforward interpretation. 
Since baryons contain valence quarks and not antiquarks, the observed excess of both protons 
and A°s over their respective antibaryons for x p > 0.2 is clear evidence for the production 
of leading baryons. The data suggest that the effect increases with x p , however more data 
are needed to study the x p dependence in detail. For x p < 0.2 the data are consistent with 
equal production of baryons and antibaryons, however the contribution from fragmentation 
is very high in this region and we cannot exclude that leading baryons are also produced at 
low x p . 

Since a meson contains one valence quark along with one valence antiquark, the inter- 
pretation of our results for mesons is more complicated. All down-type quarks are produced 
equally and with the same forward-backward asymmetry in Z° decays in the Standard Model, 
so that if a leading neutral particle such as K*° (sd) were produced equally in s and d jets 
(i.e. = — -DJf.o), then our measured D^*o would be zero. Our two highest-x p points 

are significantly positive, indicating both that there is leading K*° production and that 
more leading K*° are produced in s jets than in d jets. This is an expected consequence of 
strangeness suppression in the fragmentation process. That is, it is expected to be less likely 
for an ss to be produced from the vacuum and the s to pair up with an initial d than it is 
for a dd to be produced and the d to pair up with the initial s. 

In the case of charged hadrons such as tt~ (du), the different Z° branching ratios and 
forward-backward asymmetries of up- and down-type quarks cause a nonzero dilution of 
leading particle effects. Assuming Standard Model couplings to the Z° and equal production 
of leading ir + in w-jets and n~ in d-jets (i.e. DfL = —D™-), we calculate a dilution factor for 
our analysis cuts of 0.27. That is, we would expect to observe D n - = 0.27D^L. For purposes 
of illustration, we have fitted a line to our D p and D^o points for x p > 0.2, scaled it by the 
dilution factor 0.27, and drawn it as the dot-dashed line on figs. [He and |2~T|d. We do not 



necessarily expect that leading particle effects are identical for mesons and for baryons, but 
this line serves as a basis for a qualitative comparison. 

Our measured D n - are consistent with zero everywhere, and consistently below this line. 
This does not rule out leading pion production, but indicates that nonleading production of 
pions must be comparable or larger at all x p . This could be due to a very soft leading pion 
momentum distribution and/or a large "background" contribution of pions from decays of 
excited states such as p°, u>, rj, K*. Our measured D K - are consistently positive and above 
the line for x p > 0.2. As in the case of K*°/K*°, this indicates both production of leading 
charged kaons and more frequent production of leading K~ in s-jets than in w-jets. 
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Figure 21: Normalized differences between hadron and antihadron production in light quark 
jets. The thin dot-dashed lines in (c) and (d) represent the fit to the baryon data scaled by 
the dilution factor of 0.27 described in the text. Also shown are the predictions of the three 
fragmentation models. 
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The quantification of the total number of observed leading particles is problematic. For 
example, in the region x p > 0.2 we observe a total of 0.083±0.005 protons and 0.036±0.005 
antiprotons per light quark jet. Some of the antiprotons are expected to be "subleading" 
antiprotons produced in association with a leading baryon, since baryon number is known 
to be conserved locally [|32|, whereas others are from a non- leading baryon-antibaryon pair, 
and provide a measure of the background of nonleading protons in the high-x p sample. We 
conclude that the number of leading protons we have observed per light quark jet must lie 
between the p-p difference and the total number of protons, i.e. in the range 0.047-0.083 per 
light quark jet. Similarly, the number of observed leading A in the range 0.18 < x p < 0.5 is 
0.024-0.039. For x p > 0.26 we measure a total of 0.110±0.012 K*° and 0.023±0.010 K*° per 
light quark jet. In this case, all of these could be leading due to contributions from s and 
d jets, and so the sum gives an upper bound on the number of leading K*°/K*° produced. 
A lower bound is given by the possibility that no leading K*° are produced in d jets. In 
this case all of the observed K*° are nonleading, we expect an equal number of nonleading 
K*°, and the number of leading K*° produced is given by the K*°-K*° difference. Thus we 
have observed 0.087-0.133 leading K*°/K*° per jet with x p > 0.26. Similarly, the number 
of leading charged kaons produced in the range 0.21 < x p < 0.77 is 0.141-0.355 per jet. 

The measured normalized differences are compared with the predictions of the three 



fragmentation models in fig. 21. All models reproduce the qualitative features of our data. 
For the baryons, the HERWIG prediction drops below zero in the range in which we have 
no proton coverage; this behavior might be ruled out with more A°/A° data. The HERWIG 
and UCLA predictions rise sharply to unity at x p ~ 0.4 and are inconsistent with the proton 
data. For the mesons all models are consistent with the data. 



9 Production Ratios and Fragmentation Parameters 

Certain aspects of the fragmentation process can be studied more directly by measuring the 
relative production of two hadron species that differ by a single quantum number. We have 
calculated the ratios of differential cross sections for a number of pairs of hadron species, for 
flavor-inclusive and light-flavor events, taking into account any systematic errors common to 
the two species. The results are shown for light-flavor events in fig. |2~2]. In the cases where 



binning was different for the two hadron species in a pair, the ratio was obtained by fitting a 
curve to the denominator over a region near each x p value in the numerator. In some cases 
charged and neutral pseudoscalar kaons were averaged, and are denoted simply "K" . In all 
cases, charge-conjugate states are included in both numerator and denominator. 

The ratios of the strange mesons to pions vary rapidly with x p . In flavor-inclusive events 
(not shown), the values of each of these ratios vary over a similar range but show less 
structure, being consistent with simple powers of x p for x p > 0.04. The protompion ratio 



also varies rapidly for x p < 0.1. The other ratios shown in fig. 22 are independent of x 



■p 



within our errors. 

The K°:K ± ratio differs significantly from unity over the range 0.03 < x p < 0.09, aver- 
aging 0.86±0.03; we observe a similar difference in flavor-inclusive events (not shown), as 
has been observed previously [J]]. Assuming that primary charged and neutral kaons are 
produced equally in the fragmentation process, this implies that some hadron species is pro- 
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Figure 22: Ratios of measured differential cross sections for various pairs of hadron species 
in light-flavor events, along with the predictions of the three fragmentation models. In all 
cases the charge-conjugate states are included in both numerator and denominator. Here, 
"K" denotes the average of K°/K° and K ± . The JETSET predictions for the K*:ir+, 0:tt + , 
4>:K* and K*:K ratios have been scaled by factors of 2/3, 1/2, 4/3 and 2/3, respectively (see 
text), in order to clarify the comparison of the momentum dependence. 
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duced that decays preferentially into charged kaons. Our measured cross sections indicate 
that decays of and K* mesons would each account for only ~0.01 of the difference from 
unity. Decays of D- and 5-hadrons cannot be the source of this difference since they have 
been excluded explicitly. 

The predictions of the three fragmentation models are also shown in fig. ^2], and all de- 
scribe the qualitative features of the data. The JETSET prediction for each ratio involving 
K* or <p mesons differs from the data by a large normalization factor, and those predic- 
tions have been scaled by factors derived from fig. [15| in order to compare the momentum 
dependence with that of the data. All models underestimate the slope of the K:n + ratio, 
but reproduce those of the 0:7r + and K*:tt + ratios, overestimating the latter ratio only at 
the highest-Xp point. The x p dependence of the p:7r + ratio is reproduced by all models at 
low x p , but only by the JETSET model for x p > 0.2. However the JETSET model shows 
a normalization difference from the data of about 20%. Similar differences in the model 
predictions for the A:K ratio cannot be resolved with the current statistics. No model repro- 
duces the measured K°:K + ratio; all predict a roughly constant value of 0.98 in the range 
of our measurement. All models predict a larger value of the K*:K ratio at the highest-x p 
point than is observed in the data. A similar set of comparisons for flavor-inclusive events 
(not shown) yielded the same conclusions. 

These ratios can be used to study the suppression of baryons, vector mesons and strange 
hadrons in the fragmentation process. Quantifying such suppression at the primary frag- 
mentation level is problematic due to possible effects of different masses of the two hadron 
species in the ratio and the fact that decay products populate a different x p region than their 
primary parents. We therefore used the JETSET model, in which there are tunable param- 
eters controlling the relative production of baryons, strange hadrons and vector mesons, to 
extract suppression parameters in the context of that model. We first considered the relative 
production of pseudoscalar (P) and vector (V) mesons, traditionally expressed in terms of 
the parameter Py = V/(V + P). Since we might expect that measured ratios are not the 
same at very high x p , where leading hadron production is important, as they are lower x p , 
we defined arbitrarily a "fragmentation" region, 0.05 < x p < 0.25, and a "leading" region, 
x p > 0.45. In each region we averaged our measured K*:K ratio, and compared it with those 
obtained in the same region from the JETSET generator run with a series of input values 
of the Py parameter for strange mesons. We interpolated to find the Py value at which the 
model prediction for each ratio was equal to that measured in the data, and these values are 
listed in table [22] for the two x p regions and for both flavor-inclusive and light-flavor events. 
The two measurements in each momentum range are consistent, but the Py value measured 
in the fragmentation region is significantly higher than that measured in the leading region 
for both flavor categories. 

We next considered the relative production of baryons (B) and mesons (M), in terms of 
the parameter Pb = B/(B + M). A similar set of comparisons of our p:7r and A:K ratios 
with the predictions of the JETSET model as Pb was varied yielded the measured Pb values 
listed in table The four values extracted from the p:7r ratio are consistent. The value from 
the A:K ratio in light-flavor events is consistent with these four, but that in flavor-inclusive 
events is slightly larger. 

Information on the suppression of strangeness is available from several of our measure- 
ments. It is conventional to define a suppression factor 7 S as the probability of creating an 
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Vector:Pseudoscalar Production Parameter Py 


x p Range 


inclusive 


light-flavor 


0.055-0.219 
0.439-1.000 


0.405±0.020 
0.226±0.029 


0.433±0.033 
0.279±0.029 



Table 22: Measurements of the vector- meson fraction Py extracted from the measured 
K*:K production ratio in the context of the JETSET model. 



BaryomMeson Production Parameter Pb 


Ratio x p Range 


inclusive 


light-flavor 


p:7r ± 0.055-0.165 
A:K 0.061-0.237 
p:7r ± 0.493-0.987 


0.076±0.003 
0.101±0.003 
0.081±0.006 


0.074±0.004 
0.087±0.005 
0.081±0.009 



Table 23: Measurements of the baryon fraction Pb in the context of the JETSET model. 



ss from the vacuum, relative to that of creating a uu or dd, at a given point in the frag- 



mentation process. As has been suggested in ref. |33|, the normalized production difference 
(see section 8) at high x p between a strange hadron and its antihadron in light quark jets 
provides a robust way of investigating strangeness suppression for any neutral hadron, such 
as K*°/K*°, that is unlikely to be a decay product of a heavier primary particle. If we 
assume leading particle dominance, so that K*° can be produced only in s and d jets, and 
that the relative production in d jets is suppressed by a factor of 7 S , then we expect the 
normalized difference to be D^*o = (1 — 7 S ) /(l + 7 S ). From our point in the bin 0.5 < x p < 1 
we used this equation to derive a "direct" measurement of 7 S = 0.26 ± 0.12, where we first 
scaled our given D^,o value by 0.923 to account for the fact that we assumed contribu- 
tions from u, d and s jets in the original unfolding, whereas we now assume only d and s 
contribute. Similarly, assuming dominant production of leading K ± and accounting for the 
different branching fraction and forward-backward asymmetry of up- and down-type events, 
one expects 1.05Dx- = (1 — 0.557 s )/(l + 0.777 s ). From this we derive 7 S = 0.41 ± 0.17, 
using our D K ~ data in the range 0.47 < x p < 0.77. 

We also used the JETSET model to predict the normalized differences as a function of j s , 
and to extract from our measured _D^*o and D^- the 7 S values listed in table |24l. Also listed 



in table ^4| are 7 S values extracted in the context of the JETSET model from our measured 
K:tt + , <p:K* and A:p ratios. For each ratio, the values derived from the flavor-inclusive and 
light-flavor events are consistent. However there is a significant x p dependence in the values 
obtained from the K:tt + ratio in both flavor categories, and there are several other significant 
differences between pairs of values from the same flavor category. This indicates that the 
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Strangeness Suppression Factor, 7 S 


Ratio 


x p Range 


inclusive 


light-flavor 


D x *o 

D K - 


0.482-1.000 
0.493-0.768 


- 


0.194±0.141 
0.249±0.110 


K:ti+ 
<j>:K* 
A:p 


0.055-0.219 
0.048-0.263 
0.050-0.182 


0.236±0.016 
0.163±0.027 
0.339±0.014 


0.266±0.014 
0.184±0.052 
0.311±0.032 


K:n + 
<p:K* 


0.493-0.768 
0.482-1.000 


0.575±0.084 
0.160±0.060 


0.483±0.091 
0.239±0.075 



Table 24: Measurements of the strangeness suppression factor ^ s in the context of the 
JETSET model. The notation D h refers to the normalized differences discussed in section 
8. 



JETSET model cannot accommodate all of our data with a single 7 S value and all other 
parameters set to their default values. 

10 Summary and Conclusions 

We have measured the production of the seven hadron species 7r ± , K ± , K°/K°, K*°/K*°, 
(ft, p/p, and A°/A° as a function of scaled momentum x p over a wide range in hadronic 
Z° decays. The SLD Cherenkov Ring Imaging Detector enabled the clean and efficient 
identification of stable charged hadrons, yielding precise measurements of their production 
cross sections, as well as the identification of relatively clean samples of the strange mesons 
K*°/K*° and <fi reconstructed in decay modes containing charged kaons. Our measurements 
of differential production cross sections, total cross sections and ratios of production of these 
hadron species in flavor-inclusive hadronic Z° decays are consistent with averages of those 
from experiments at LEP. 

Using the SLD vertex detector to isolate high-purity light- and 6-tagged event samples, we 
have measured the production of these seven hadron species in light-, c- and 6-flavor events. 
Significant differences between flavors were found, consistent with expectations based on the 
known properties of B and D hadron production and decay. Our 7r ± , K ± and p/p data 
at high x p were used to test the predictions of Gribov and Lipatov for the shape of the x p 
distribution of primary leading hadrons as i p 1. We find the predictions of the theory to 
be consistent with the flavor-inclusive (light-flavor) meson data for x p > 0.66 (x p > 0.47) and 
with the proton data for x p > 0.43 (x p > 0.38). The shape of the £ = — ln(x p ) distribution for 
each hadron species in events of each flavor is consistent with the Gaussian form predicted 
by MLLA QCD+LPHD near its peak. The peak positions for each hadron species in 
light-flavor events are more consistent with a monotonic dependence on hadron mass than 
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those in flavor-inclusive events. 

Using the large forward-backward asymmetry induced by the polarized SLC electron 
beam to separate light quark from light antiquark hemispheres, we have updated our mea- 
surements of hadron and antihadron production in light quark jets. Differences are observed 
at high x p between baryon and antibaryon production, which is evidence for the production 
of leading baryons, i.e. baryons that carry the quantum numbers of the initial quark. Dif- 
ferences are also observed for both pseudoscalar and vector i^-mesons, which indicate not 
only leading production of these two hadron species but also that leading strange mesons 
are produced more often from initial s quarks than from initial u or d quarks. 

Our data were used to test the predictions of three fragmentation models with default 
parameters. In most cases these simulations reproduced the data to within a few percent. 
However the JETSET 7.4 model predicts too many p/p, K*°/K*° and 4> mesons at all x p , and 
too many K and K°/K° at low x p . The UCLA model predicts too many pions in the 2-20 
GeV/c range, a shoulder in the x p distributions for baryons at high x p , and larger differences 
between baryon and antibaryon production at high x p than are seen in our light-quark data. 
The HERWIG 5.8 model predicts a shoulder in the x p distribution for most hadron species 
at high x p , a large excess of low-x p pions and kaons in b- flavor events and of medium-Xp pions 
in c-flavor events, and a rapid variation in the baryon-antibaryon differences as a function of 
x p . All models predict a charged:neutral kaon ratio very close to unity, which is inconsistent 
with our light-flavor and flavor-inclusive data. Also, no model is consistent with the x p 
dependence of either our K:ir ratio or our K*:K ratio. 

We have studied several parameters of the fragmentation process. The differences be- 
tween kaon and antikaon production in light quark jets allow two new, direct measurements 
of strangeness suppression at high momentum. We have also used our ratios of production of 
pairs of hadron species to extract fragmentation parameters in the context of the JETSET 
model. We find the vector:pseudoscalar meson parameter to be dependent on x p , and the 
strangeness suppression parameter to be dependent both on x p and on the hadron species 
used to form the ratio. 
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